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ORBITAL RECONNAISSANCE PROGRAM 
FOR VIKING CLASS MISSIONS 


By Richard N. Green and Sue W. Souders 
Langley Research Center 

SUMMARY 

A computer program has been developed to meet the needs of Project Viking in the 
area of orbital reconnaissance. This program is presented as an analytical tool for pre- 
liminary mission design and analysis. 

The program has the capability of investigating the time history of both the space- 
craft groundtrack and the camera footprint on the rotating surface of Mars. Various 
options are available to control where and how the photographs are taken. A sequence of 
photographs can be taken over a desirable lighting band, between two orbital true anoma- 
lies, or between latitudes of interest. Within the interval of desirable photography the 
photographs can be taken on either a time increment or a true anomaly increment. Also 
included is the option to take photographs on an overlap area consideration. Both vertical 
and nonvertical photography are available as program options. A description of the pro- 
gram input and output, a FORTRAN listing of the program, and samples of the input and 
output are included. 


INTRODUCTION 

A computer program, PANDG (Photography and Groundtrack), has been developed to 
meet the needs of Project Viking in the area of orbital reconnaissance. Specifically, 
PANDG was developed in support of the Langley Viking Project Office for use in prelimi- 
nary mission design and analysis. 

The program has the capability of investigating the time history of both the space- 
craft groundtrack and the camera footprint on the rotating surface of Mars. Various 
options are available to control where and how the photographs are taken. A seqi'ce of 
photographs can be taken over a desirable lighting band, between two orbital true >- 
lies, or between latitudes of interest. Within the interval of desirable photography the 
photographs can be taken on either a time increment or a true anomaly increment. Also 
included is the option to take photographs on an overlap area consideration. Both vertical 
and nonvertical photography are available as program options. 



The program has been designed to output either photographic data or groundtrack 
data or the two in combination. The groundtrack output consists of the Julian date, the 
calendar date, and the time elapsed from the initial true anomaly. Also included are the 
latitude and longitude of the subsatellite point, the altitude, inertial velocity, true anomaly, 
and flight -path angle of the spacecraft. Additional parameters include the lighting angle 
at the subsatellite point, the hour angle of the Martian vernal equinox measured from the 
zero meridian, the argument of periapsis, and the longitude of the ascending node of the 
spacecraft orbit. Equations for most of these parameters are developed herein. The 
photographic output consists of parameters of interest at the four corners of the camera 
footprint, at the midpoint of each of the four sides, and at the center point of the photo- 
graph. At each of these nine points the following parameters are output: latitude, longi- 
tude, lighting angle, slant range distance, and static resolution. One of the most important 
photo output parameters is the ratio of the overlap area between the present camera foot- 
print and the previous footprint. The development of these parameters is also presented. 

One advantage of PANDG is the capability to control the overlap area. Previously, 
some investigators have used forward overlap and side overlap as a measure of the over- 
lap area; this definition was commonly used for Lunar Orbiter photography. However, 
for Mars photography, a more complete definition is needed because of the geometry 
created by a rapidly spinning planet. Therefore, overlap was defined as the ratio of the 
overlap area to the area of the previous photograph. For simplicity, the sides of the foot- 
print were considered great circles instead of minor circles. An option is available to 
control the overlap between consecutive photographs. 

The information necessary for the implementation of the program is contained in 
the appendixes of this paper. The function of the main program is outlined with a flow 
diagram; the purpose of each subroutine is set forth in appendix A. A FORTRAN listing 
is included in appendix B, in addition to sample input and output in appendix C. 

SYMBOLS 

X vector from center of planet to corners of previous photographic footprint 

on surface, kilometers 

a semimajor axis, kilometers 

B vector from center of planet to corners of present photographic footprint 

on surface, kilometers 

C component of unit C 


2 


c 

c 

D 

D 

d 

E 

e a 

E C 

e 

F 

f 

H 

H a 

H p 

h x ,hy,h z 

h 

t 

i 

JD 

J 20 


camera vector directed from spacecraft along one corner of photograph to 
surface, kilometers 

unit vector from spacecraft along one corner of photograph 

component of D, kilometers 

arbitrary vector within overlap area, kilometers 

magnitude of C, kilometers 

spherical excess, degrees 

spherical excess of footprint formed by A^,A2,A3,A4 vectors, degrees 

spherical excess of overlap area, degrees 

eccentricity 

vector from center of planet to one corner of overlap area, kilometers 

true anomaly, degrees 

altitude, kilometers 

altitude of apoapsis, kilometers 

altitude of periapsis, kilometers 

components of h 

unit vector in azimuth direction (heading) 

footprint intersection vector, kilometers 

inclination, degrees 

Julian date, days 

second zonal harmonic of planet 


n x ,n y» n z 


components of n 



»-s> 


n 

p,Q,w 


r 


r s 


r x> r y> r z 


r 

S 

S 


t 

V 


V/H 

X,Y,Z 

x' ,y',z' 

x,y,z 
x’ ,y' ,z' 


x 

4 


unit vector normal to orbital plane 

P points toward periapsis, Q is in orbital plane advanced to P by a 
right angle in direction of increasing true anomaly, and W completes 
right-handed system 

magnitude of r , kilometers 

radius of surface, kilometers 

components of r 

unit vector in radial direction 

radius vector, kilometers 

component of unit vector directed from planet toward Sun 

unit vector from center of planet parallel to arrival asymptote of 
incoming hyperbola 

time, days 

velocity, kilometers/second 

hyperbolic excess velocity in areocentric coordinate system, 
kilometers/ second 

horizontal velocity relative to surface divided by altitude, second" * 
Cartesian coordinate axis system 
primed Cartesian coordinate axis system 
rectangular Cartesian coordinates, kilometers 

rectangular Cartesian coordinates in the primed system, kilometers 
(see sketch (h)) 

spacecraft position vector, kilometers 


a 


longitude, degrees 


0 angle between spacecraft position vector and camera vector, degrees 

(see sketch (e)) 

0* critical camera angle, degrees (see sketch (e)) 

y flight -path angle, degrees 

AF magnitude of AF 

AF = F - D 

6 latitude, degrees 

rj angle between sides of overlap area, degrees (see sketch (i)) 

8 hour angle, degrees 

© longitude of D, degrees 

X angle between radius and camera center line, degrees (see sketch (b)) 

Id gravitational constant, kilometers^/second^ 

£ angle between sides of camera footprint, degrees (see sketch (f)) 

S azimuth, degrees 

$ latitude of D, degrees 

0 angle between radius vector and Sun vector, degrees 

O longitude of ascending node, degrees 

co argument of periapsis, degrees 

camera forward half-angle, degrees (see sketch (b)) 
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camera side half-angle, degrees (see sketch (b)) 


^s 

Subscripts: 


o initial conditions 

R measured in rotating coordinate system 

r,h,n rectangular Cartesian components in r,h,n coordinate system 

x,y,z rectangular Cartesian components in X,Y,Z coordinate system 

x',y' ,z' rectangular Cartesian components in primed system of coordinates 
(see sketch (h)) 


Caret (~) over a symbol denotes unit vector. 

Dot over a symbol denotes derivative with respect to time. 

Numerical subscripts are used to distinguish similar quantities. (For example, 
A^,A2,A3,A4 are the four vectors defined as A.) 

METHOD OF CALCULATION 


The PANDG computer program has been designed to compute areographic photo- 
graphic and groundtrack data rapidly. The motion of the spacecraft about Mars has been 
considered Keplerian motion with the inclusion of the second zonal harmonic J20, which 
perturbs only the longitude of the ascending node and the argument of periapsis. The 
program is complicated only by the variety of options which are available to the user. 
These options, however, are necessary to provide the flexibility needed for mission anal- 
ysis. The program operates in three modes. The first mode produces only groundtrack 
data, the second produces only photographic data, and the third mode combines the two 
types of data. With the third option the groundtrack of a number of revolutions will be 
output with photographic output during regions of interest. 

The initial state of the spacecraft can be input in one of three ways. The most 
obvious definition of the initial state consists of the six Keplerian orbital elements a, e, 
i, u>, ft, and f. Often, however, the altitudes of periapsis and apoapsis are more conve- 
nient to input. Therefore, the option to replace a and e with H a and Hp has been 
included. The third option derives the initial state from a knowledge of the incoming 
approach hyperbola. (See ref. 1 .) To exercise this option, the right ascension and 
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declination of the approach asymptote S must be input along with the hyperbolic excess 
velocity and the desired inclination of the spacecraft orbit. In addition, the shape 

of the orbit must be defined by H a and Hp. This option assumes a deboost from a 
hyperbolic periapsis to an elliptical periapsis and defines the initial true anomaly of the 
spacecraft as zero. Since PANDG calculates data relative to the rotating surface, the 
initial time is important. Allowance has been made to input either an initial Julian date 
or an initial calendar date. 

From the initial state, the spacecraft is stepped either to the point of the first 
groundtrack output or to the point of the first photograph, depending on which occurs first. 
Many times it is desirable to increment the spacecraft in time. Often, however, a true 
anomaly increment seems more appropriate. Therefore, both types of increments have 
been included as options. 

If the program is exercised in one of the two photographic modes, the region over 
which photographs are to be taken must be specified. This region can be a Sun angle 
band, the section along the trajectory defined by two true anomalies, or a latitude band. 
Also included is the option to take a single photograph at a given latitude. Within the 
photographic region of interest, multiple photographs are taken. The spacing of these 
photographs can be specified in one of two ways, that is, either on the groundtrack step 
(true anomaly or time) or on an overlap area consideration. In addition to the above 
specifications, the photography must be defined as either vertical or nonvertical. If the 
nonvertical option is employed, the camera axis is rotated about the heading direction 
from a vertical direction through an angle X. In other words, the camera axis is rotated 
out of the orbital plane by an amount of A. 

Termination of a specific case can be defined by three stop options - time, orbit 
number, or the central angle traversed by the spacecraft. There is no limit to the period 
of time which can be investigated. As the spacecraft position is incremented in time, the 
position of Mars is also updated with the mean orbital elements of the planet about the Sun. 
The location of the Martian prime meridian is also updated so that areographic param- 
eters, expressed in the Mars centered rotating coordinate system, can be examined. 
Therefore, the number of days that are investigated is at the discretion of the user. 

Initial State of the Spacecraft 

Dependent on which input option is exercised, the six Keplerian orbital elements of 
the spacecraft are calculated. The first option calls for all six elements. The second 
option expresses the shape of the orbit in terms of H a and Hp, the altitudes of apo- 
apsis and periapsis, respectively. For the second option, the semimajor axis a and 
the eccentricity e are given by 
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H a + Hp + 2r s 
2 


a 

where r g is the radius of Mars. The third option derives the orbital elements from 
the arrival hyperbola. (See ref. 1.) The six orbital elements are then converted to 
inertial Cartesian components. 


Groundtrack Parameters 

The groundtrack output consists of a number of parameters, some of which relate 
to the Martian surface. To calculate these parameters, the state of the spacecraft must 
be expressed in the rotating equatorial coordinate system. The relationship between the 
two Mars centered coordinate systems, the inertial equatorial coordinate system (areo- 
centric) and the rotating equatorial system (areographic), is through the hour angle of the 
Martian vernal equinox as measured from the prime meridian. The hour angle 6 is 
given by (ref. 2) 

e = e 0 + e At (l) 

where 

e Q = 149°. 475 
6 = 350.891962 deg/day 
At = JD - 2418322.0 

The present time is expressed in terms of Julian date JD. 


I 


! 


Z 



Sketch (a) 
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It is now an easy task to express the spacecraft state in the areographic coordinate 
system. From sketch (a) and reference 3, 

~\ 

x R = x cos 9 + y sin 6 
y R = -x sin 0 + y cos 6 


Z R ~ z 


x R = x cos 9 + y sin 9 + y R 0 
y R = -x sin 9 + y cos 6 - x R 0 




Z R = z 


The longitude of the subsatellite point is given by 

yR 






sm a = 


cos ot - 


^R 2 + yR 2 
X R 




yxj-2 v_2 


{ r + yR 

where 0° £ a £ 360°, and the latitude is given by 

‘- 1 ® 

where -90° = 5 s 90° and 


( 2 ) 


(3) 


(4) 


The altitude is simply 


-i 


R 


+ y R 2 * 


J R 


H = r - r s 

and the inertial velocity of the spacecraft is given by 



To find the inertial flight -path angle y, the vertical component of the inertial veloc- 
ity vector, which is found by the dot product, is divided by the velocity, that is 


y = sin 


-i f xx + yy + zz \ 
V rv J 


where -90° = y = 90°. The Sun angle at the subsatellite point is defined as the angle 
between the spacecraft radius vector and the vector to the Sun. When the Sun is directly 
overhead, the Sun angle is zero. To obtain the unit Sun vector, a series of subroutines 
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which utilize the mean orbital elements of Mars about the Sun is exercised. A discus- 
sion of this procedure is given in reference 1. Once the Sun vector has been obtained, 
the vector dot product yields the Sun angle at the subsatellite point, that is 

,/x S + yS + zS \ 

<p = cos -1 ! — — -2 -J (5) 

where 0° s 0 % 180° and (S x ,S y ,S z ) is the unit vector directed toward the Sun. Also 
included in the set of groundtrack parameters is V/H, or the horizontal velocity relative 
to the surface divided by the altitude of the spacecraft. The horizontal velocity can be 
found by considering the magnitude of the vector cross product of the relative velocity 
vector and the radius vector, that is 


V(horizontal) = 



= V R sin(90° - y) 


where (90° - y) is the angle between the two vectors, or the complement of the relative 
flight-path angle. Therefore, V/H is given by 


V/H = 


_ [( y R Z R ~ ^R Z r) 2 + ( X R Z R ~ x R & r) 2 * ( x R y R ~ x R y It) 2 _ 


1/2 


rH 


In addition to the above groundtrack parameters are the following: Julian date, calendar 
date, orbit time, orbit number, and the true anomaly of the spacecraft. 


Photographic Parameters 

The main objects of interest at a photographic point are the camera footprint on the 
surface (that area of the surface photographed by the onboard camera) and the set of 
parameters which are evaluated at each of the four corners of the footprint and at the mid- 
point of each side. The location of the footprint, however, depends on the camera system 
and its characteristics. The camera system discussed herein is defined as being alined 
with the inertial orbital plane and is described by a forward half -angle and a side 
half-angle (See sketch (b).) In addition, the capability to rotate the camera center 

line out of the plane of motion by an angular amount X has also been included. For verti- 
cal photography, the camera center line is directed downward along the radius vector 
(X = 0). To facilitate the calculation of the points where the four corner vectors pierce 
the surface, consider the spacecraft-centered coordinate system r,h,n (sketch (c)). The 

A 

unit vector ? is directed along the radius vector to the spacecraft, and the vector h 
(heading) is in the azimuth direction. The unit vector n is normal to the plane of motion, 
that is, n = r x ft. The unit vector C directed from the spacecraft along one of the cor- 
ners of the photograph, defined by a positive forward half-angle i/'f and a positive side 
half-angle i// s , is illustrated in sketch (d). Also, a positive out-of -plane angle X has 
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Sketch (h) 



Sketch (c) 


n 



Sketch (d) 
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been considered. The components of the corner vector in the r ,fi,n system are there- 
fore given by 

C r = -cos (X + xjy s ) cos i//j| 

C h = cos (X + \p s ) sin i// f ) (6) 

C n = sin (X + i/'s) 

The corner vector in the X,Y,Z system remains to be expressed. The unit vec- 
tor r expressed in the X,Y,Z system is 


f = f - ( r x’ r y> r z) 


and the unit vector n is given by 


- rxx _ /„ „ _\ 
n — 7 — tt - ( n x» n y> n z) 
r x x 


The unit vector h is 

h = n x r = (h x ,hy,h z ^ 
Therefore, the unit vector C in the X,Y,Z system is 


C z 


n s 

n. 


r y h y 


r z h z n z 



Cr 


Ch 


_ c n 


From sketch (e) , the angle /3 between the radius vector x and the corner vector C 
is given by 


j3 = cos 


-1 


(-x . c) 


where 0° = 0 s 180°. If 0 exceeds the value of 0*, the corner vector does not pierce 
the surface. From sketch (e), 

t - sm -*($ 

where 0° = 0* i 180°. Hence, for the corner vector to pierce the surface, 0 < 0*. The 
slant range distance d is given by the law of cosines as 

r s 2 = r 2 + d 2 - 2rd cos 0 
or 


2r cos 0 • ^4r 2 cos 2 0 - 4^r 2 - r g 2 j 
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Since the extended corner vector pierces the surface at two points and since the smaller 
of the two solutions for d is the desired solution, the radical is negative and 


d = r cos /3 - yr^cos^ - r^ + r g ^ 

Therefore, the point at which the camera corner vector pierces the surface of the planet 
is given by 

F = x + C 


where 

C = dC 


The other three corners are found in a similar manner by considering ±\p s and ±\p^ 
in equation (6). The latitude and longitude of the corner vector F are found with equa- 
tions (1) to (4) and the Sun angle with equation (5). 


Camera Footprint Overlap 

Another photographic parameter of interest is the camera footprint overlap. The 
camera overlap area, as defined in PANDG, is the surface area contained in both the pres- 
ent camera footprint and the previous camera footprint. The ratio of this overlap area to 
the area of the previous footprint is designated as the camera footprint overlap or simply 
as the overlap. This parameter is important in designing a sequence of photographs so 
that they can be overlaid to form a map of the area photographed. In addition, overlap is 
necessary to obtain stereo effects from the photographs. Previously, some investigators 
have used forward overlap and side overlap as a measure of the overlap area. This 


r 
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definition was commonly used for Lunar Orbiter photography. However, for Mars photo- 
graphy, a more complete definition is needed because of the geometry created by the 
higher spin rate of the planet. In certain instances the surface with the camera footprint 
rotates faster than the spacecraft revolves in its orbit. Therefore, the footprint is ahead 
of the spacecraft, and forward overlap would not be adequate. 

The camera footprint is that area of the surface photographed by the onboard cam- 
era. As shown previously, the four corner points or vectors can be found. The sides of 
the footprint are actually arcs of minor circles, which are difficult to analyze mathemat- 
ically. To simplify computations, the sides have been considered as great circles. This 
approximation is justifiable for photography taken at low altitudes (small footprint). 
However, as the footprint becomes large, the error associated with this approximation 
becomes more pronounced. Since most of the photography will be at low altitudes, this 
simplification seems reasonable. 

The overlap is defined by eight corner vectors. (See sketch (f).) The four vectors 
Aj, A2, A3, and A4 represent the four corners of the previous footprint. Similarly, 
the present footprint is defined by the four B vectors. For the particular geometry 
shown, the overlap area has four corners which are typed as either interior vectors or 
intersection vectors. Corner vectors of the overlap area which are also corner vectors 
of one of the footprints are classed as interior vectors. Examples of these are Bj and 
A3. An intersection vector occurs at the corner formed by the intersection of one side 
from each of the two footprints. Examples of these are Tj and I2. 

The surface area of the previous footprint is found by considering the number of 
spherical degrees in the footprint. Since the entire spherical surface of Mars contains 
720 spherical degrees and since the surface area is easily obtainable, the area of the 
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footprint follows from the number of spherical degrees it contains. However, to obtain 
the overlap, which is a ratio of areas, only the computation of the spherical degrees con- 
tained in the area is necessary and not the area itself. The number of spherical degrees 
in an area is equal to its spherical excess. By definition, the spherical excess of a spher- 
ical polygon is the difference between the sum of its angles and the sum of the angles of a 
plane polygon having the same number of sides. (See ref. 4.) Therefore, the spherical 
excess of the footprint formed by the A vectors is given by 

= (£j + ^2 + ^3 + ^ 4 ) " 

The ij angles are the angles between the great circular arcs which connect the corner 
points. These angles are most easily found by considering the vector normal to each 
side and computing the angles between the normals. Therefore, is given by 


= cos' 


x ^4) * x ^2) 


Aj XA 4 Aj x A 2 | 


where 0° s < 180°. The other three i; angles are found in a similar manner. 

Prior to the calculation of the spherical excess of the overlap area, the corner vec- 
tors must be found. First, the B vectors are examined to determine if they are interior 
vectors. If Bj is an interior vector, the sum 
of the spherical excesses of the four triangles 
Bi-Ai-A 2 , Bj-A2-Ag, BJ-A3-A4 , Bj-A 4 -Aj 
will be equal to the spherical excess of the foot- 
print formed by A vectors. However, if these 
excesses are not equal, the vector considered is 
not an interior vector but lies outside of the foot- 
print area. (See sketch (g).) The spherical 
excess of the triangle B 4 -Aj-A 2 is given by 



E = 


cos 


-1 


(^1 x %) • (A t x A 2 ) 
| X 1 x5 l|| X l xX 2| 


1 n ac “ 1 

(a 2 x A x ) • (A 2 X Bi) 

j. PnQ ~ 1 

(5l X A 2 ) . (bj x A i) 

T CUD 

|a 2 x Aj| 1 A 2 x Bj| 

T vUQ 

|bi xA 2 ||b! xA t | 


The spherical excesses of the other three spherical triangles are found in a similar man- 
ner, and the sum indicates whether Bj is an interior point. This sequence is then 
repeated for both the A and the B vectors until all interior vectors are found. 
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Once the interior vectors are found, the geometry is examined to determine if any 
intersection vectors exist. First, the side joining Aj and A2 is examined to deter- 
mine if any of the sides formed by B vectors intersect it. The intersection vector is 
found by the vector cross product of the normals, that is 

T = ±(Aj x A 2 ) x (gj x B 4 ) 

where T is the intersection vector and is multivalued since great circles intersect at 
two points. For T to be a corner vector of the overlap area, it must lie between A\ 
and A 2 and also between Bi and B4. If the sum of the angles between Aj and I 
and between I and A 2 equals the angle between Aj and A 2 , then I is surely 
between Aj and X 2 * The same must be true for the B vectors. Therefore, for I 
to be an intersection vector, 


and 


cos 


-1 


A x • 

T 

_l pnc " 1 

T ■ a 2 


A 1 * a 2 


T| 

T UUo 

T a 2 

— UJD 

mn 


cos 


Bi . I 


Bi 


+ cos 


I • B/ 


b 4 


= cos 


-1 


Bl 

•b 4 

Bl| 

|b 4 [ 


Both values of I must be examined as candidates for an intersection vector. As can be 
seen in sketch (f) , I satisfies the first of the two criteria but fails the second. There- 
fore, I is not a corner vector of the overlap area. In turn, each of the four sides of the 
footprint formed by B vectors is examined in combination with the four sides of the foot- 
print formed by A vectors. For the geometry of sketch (f) , the intersection vectors Tj 
and I 2 are found in this manner. 


Once all the corner vectors of the overlap area are found, the overlap area is deter- 
mined. However, the interior and intersection vectors found are not necessarily in clock- 
wise or counterclockwise order. To facilitate the calculation of the overlap area, the 
corner vectors are ordered by considering the azimuth angle of each vector referenced to 
an arbitrary vector within the overlap area. The corner vectors are redefined as Fj 
where i = 1, 2, 3, . . ., n, and the reference vector D is formed by averaging the first 
two corner vectors, that is 


The vectors from D to each of the corner vectors are given by 

AFj = F A - D 


where i = 1, 2, 3, . . ., n. 
16 



The azimuth angle is found by considering the geometry of sketch (h). The primed 
axis system is alined such that the X'-axis is in the radial direction, the Z'-axis is per- 
pendicular to X’ and in a northerly direction, and the Y'-axis completes the triad. The 
rotation from the unprimed system to the primed system is accomplished by 


1 

X 

= 

cos $ 0 sin $ 


cos © sin © 0 

X 

T 

y 

= 

0 10 


-sin © cos 0 0 

y 

t 

Z 

= 

-sin $ 0 cos $ 


0 0 1 

z 


where 


sin 0 = 


cos © = 


sin $ = 




+ D Z 


2 


cos $ = - sin^$ 
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Therefore, AF can be expressed in the primed system as 


AF X , =AF X cos <f> cos 0 + AFy cos 3> sin 0 + AF Z sin 
AFyt = -AF X sin 0 + AFy cos 0 

AF z t = -AF X sin $ cos 0 - AFy sin <f> sin 0 + AF Z cos <i> 


and the azimuth angle 2 is given by 

sin 2 = 

cos 2 = 


AFy, 


]/a Fy , 2 + AF z , 2 

AF z’ 
^AFy , 2 + AF Z , 2 


where 0° ^ 2 % 360°. In this manner the azimuth angle for each of the corner vectors 
AF i is found and ordered such that they are clockwise (sketch (i)). 



Sketch (i) 


The spherical excess Eq of the overlap area is now given as 

E C = (^i + V 2 + *? 3 + *? 4 ) - 360° 

where the rj angles are found in the same manner as the £ angles. The overlap, or 
the ratio of the overlap area to the area of the previous footprint, is simply 

E r 

Overlap = — — 

E A 

where the range of the overlap is from 0 to 1. 
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The number of corner vectors comprising the overlap area is not necessarily four 
but can vary from three to eight. Examples of various geometries are shown in sketch (j). 
Although the method described was derived by considering a four-cornered overlap area, 
which is the most frequent case, the procedure is applicable for all geometries. The 
computer program, PANDG, has the capability to calculate the overlap for all six pos- 
sible geometries. 




Four corners 



Five corners 





Sketch (j) 


RESULTS AND DISCUSSION 

The equations developed herein and the program options outlined have been incor- 
porated into the computer program designated PANDG. The program was written in 
FORTRAN IV language for a digital computer and contains a main program and 31 sub- 
routines. The information necessary for the implementation of the program is contained 
in the appendixes of this paper. Figure 1 is an example of the type of plots which can 
result from the execution of PANDG. 
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Figure 1.- Photography footprints between latitudes of -25° and 30° 
and groundtrack for a satellite orbit about Mars. 

The groundtrack of the spacecraft has been traced by the subsatellite points at 10° 
intervals of true anomaly. Also shown is a series of camera footprints which resulted 
from the choice of several program options. The region of photographic coverage was 
designated as being between the latitudes of -25° and 30°. Within this region, vertical 
photography was taken on an overlap consideration. The first in the sequence of foot- 
prints occurred at the point where the camera center line crossed the 30° latitude line. 
Photographs were then taken such that a 25-percent overlap existed between consecutive 
frames. Photographic coverage continued until all four corners of the footprint were out 
of the latitude band. PANDG then continued to output groundtrack data until one of the 
program stop conditions was reached. For the particular case considered herein, the 
program stop occurred when 400° of central angle had been traversed by the spacecraft. 
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CONCLUDING REMARKS 


A computer program has been developed to meet the needs of Project Viking in the 
area of orbital reconnaissance. The program is intended as an analytical tool for pre- 
liminary mission design and analysis. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 2, 1970. 
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PROGRAM DESCRIPTION 


The program PANDG has been written entirely in FORTRAN IV computer language 
for the Control Data 6600 digital computer and contains a main program and 31 subrou- 
tines. The modular form in which the program was written resulted in a field length 
(storage requirement) of 60 OOOg. A brief statement of the purpose of each subroutine 
contained in the program follows. 


CALJUL 

CARCON 

CARSPH 

CONCAR 

DOT 

E EARTH 

EMARS 

EULER 

GPRINT 

JULCAL 

LATLNG 

LOGIC 

OPRINT 

ORBIT 


Converts calendar date to Julian date 
Converts Cartesian coordinates to conic elements 
Converts Cartesian coordinates to spherical coordinates 
Converts conic elements to Cartesian coordinates 
Calculates the angle between two vectors 

Calculates the mean heliocentric position and velocity of Earth 

Calculates the mean heliocentric position and velocity of Mars 

Performs an Euler rotation 

Calculates and writes groundtrack data 

Converts Julian date to calendar date 

Converts Cartesian position to latitude and longitude 

Examines program options to determine whether the next event is a 
groundtrack or photography event and calculates the time increment 
and true anomaly increment to the next event 

Writes program inputs 

Calculates Keplerian orbital elements for a periapsis-to-periapsis deboost 
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OVERLAP 

PPRINT 

PRATIO 

PRECES 

QADRAT 


Calculates the time increment to the next camera footprint for a given 
overlap ratio 

Calculates and writes photography data 
Calculates the overlap ratio between two footprints 

Transforms mean Earth equinox and equator coordinates from one epoch 
to another epoch 

Solves an equation of the form AX 2 + BX + C = 0 for the real roots 


RAY 

RECEQ 

REQMEQ 

REQPEQ 

RXYZPQW 

STDCASE 

SUNBAND 

TCONIC 

TINVS 

UPDATE 


Calculates the point at which the camera vector pierces the surface 
of the planet 

Rotates a vector from the coordinate system of the mean equinox and 
ecliptic of date to the coordinate system of the mean Earth equinox 
and equator of date 

Rotates a vector from the coordinate system of the mean Earth equinox 
and equator of date to the coordinate system of the mean Mars equinox 
and equator of date 

Rotates a vector from the coordinate system of the mean Earth equinox 
and equator of date to the coordinate system of the mean planet equinox 
and equator of date 

Rotates a vector from the X,Y,Z to the PQW coordinate system 

Defines program constants and standard input data 

Calculates the two positions in orbit which correspond to a given lighting 
angle 

Calculates the time from periapsis passage for a given true anomaly 

Converts mean anomaly to eccentric and true anomaly 

Updates all time dependent parameters such as true anomaly, hour angle, 
and so forth 
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VECTOR Calculates the positions of the Sun, Earth, and Canopus in the coordinate 
system of the center of the planet and the planet equator 

WELLS Updates the argument of periapsis and the longitude of the ascending node 

The main program is relatively simple, since its main purpose is administration. 
The subroutine STDCASE is called first to define the program constants, to initialize 
various parameters, and to define any standard data which it might contain. Afterwards, 
the input data are read into the program by means of a FORTRAN IV namelist. Only two 
options are considered in the main program. The input option (IPUT = 1, 2, 3) denotes 
the set of parameters input to describe the initial state of the spacecraft. The other 
option (IDATE = 1,2) denotes which date was used; the initial date can be input as either 
a calendar date or a Julian date. Dependent on the options, the initial orbital elements 
and the initial Julian date are calculated. In addition, the initial hour angle HA and 
the unit Sun vector (subroutine VECTOR) are calculated. Afterwards, subroutine OPRINT 
is called to write all the program inputs. 

At this point the program is ready to perform its primary task of calculating photo- 
graphic and groundtrack data. Subroutine LOGIC is called to determine whether the next 
event is a groundtrack event or a photographic event. This information is contained in 
the output parameter JJ. The time increment TSTEP and the true anomaly increment 
FSTEP to the next event are also determined inside LOGIC. Next, subroutine UPDATE 
is exercised to update all time-dependent parameters according to TSTEP and FSTEP. 
Before processing the next event, however, the program stop conditions are examined. 

If the time has exceeded the time stop or the central angle has exceeded its limit or the 
orbit-number stop has been reached, then the case is terminated and a new set of input 
data are read. If none of these conditions are met, then either the groundtrack data or 
the photographic data are calculated and output. The program then proceeds to determine 
the next event. This cycle is repeated until the case is complete. A flow diagram of the 
main program follows. 
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Call STDCASE 
t = 0 
f = 0 
ONtM = 1 


Read input data 


a = (H p + Ha + 2r s )/ 2 

e = 1 - ( r s + 


A 

r ipot 
\ / 

^ — 

A 

X 

1 



Call ORBIT 


Calculate angular 
rate and period 


Call CALJUL 



Call JULCAL 


HA = HAREF + HADOT (JD-JDREF) 


Call VECTOR 


Call OPRINT 


Initialize parameters inside sub- 
routine LOGIC 

Call LOGIC (0, TSTEP, FSTEP, JJ) 
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The program logic is contained in the subroutine LOGIC. The purpose of the sub- 
routine is to examine the program options in order to determine whether the next event 
is a photographic event or a groundtrack event. In addition, it calculates the time and 
true anomaly increments to the next event. The logic involved in this calculation is 
rather complicated, as can be seen from the flow diagram of the subroutine. A number 
of different options are examined in LOGIC. The option IPG determines whether the 
program generates groundtrack data (IPG = 1) , photographic data (IPG = 2) , or the two 
in combination (IPG = 3). If groundtrack data are generated, then ISTEP determines 
whether this data is output on time (ISTEP = 1) or true anomaly (ISTEP = 2). If photo- 
graphic data are generated, the camera system is defined as either vertical (IPHOTO = 0) 
or nonvertical (IPHOTO = 1). For a single photograph, ISINGLE = 1. Usually, however, 
a series of photographs is taken. The region of photographic interest is defined by the 
I FRAME option. Within this region the photographs can be taken on either an overlap 
consideration (10 LAP = 1) or a constant step (10 LAP = 0) . Subroutine LOGIC considers 
all these different options to determine the next event. The flow diagram follows. 

Flow Diagram for 

Subroutine LOGIC 
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PROGRAM LISTING 

Program PANDG contains a main program and 31 subroutines. A FORTRAN IV 
listing of the program is contained in this appendix. The purpose of each subroutine is 
stated in the listing. Several of the subroutines contained in PANDG are not presented 
herein, since their listings appear in reference 1. The following subroutine listings 
have beem omitted: CARSPH, CONCAR, DOT, EEARTH, EMARS, EULER, JULCAL, 
LATLNG, ORBIT, PRECES, QADRAT, RECEQ, REQMEQ, REQPEQ, RXYZPQW, 
SUNBAND, TCONIC, TINVS, VECTOR. 
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PROGRAM PANDG t INPUT , OUTPUT, TfPE5 = INPUT.TAPE6=0UTPUT> 

C. 

COMMON/BLCCK1/ JDC, JD.DATE,CNUM,»-NOL AST , DEL T , DEL F 
CQMMON/BL QCK2/ A * E*XI»W»C»F,PERICO»ANGRAT»FO*FTOTAL»TTOTAL 
COMMON/RLOCK3/RD,CR.U,RS,P I ,H AR EF, HADOT , JDREF 
COMMON /BL0CK4 / SX * SY» S Z , EX ,EY »E Z *CX*CY»CZ 

COMMCN/BLOCK5/IPG, IPUT, I CATE, I STEP, IOCC , I PHOTO , I MANU , IB , I OL AP , I FR A 
♦ME.ISINGLE 

COMMON/ BL 0CK6/T STOP *FSTCP,CSTOP 

COMMON/ BL0CK7/H AC ,HPO , AC,EC,XIC,WC» 00, V INF, SDEC,SRA, BETA 
COMMON/BL 0CK8/XJ 20, RES, PH II, PHI 2, CVLAP.CANGF *C ANGS *F PHOT 01 ,FPHCTC2 
*.FOOTL,L AMPC A,L ATI , L AT2 » PL AT 
DIMENSION D ATE ( 6 > »F00TL<3»4) 

REAL LAT1,LAT2,L AMBDA, L ATL AST 
REAL JD,JD0, JDREF 
INTEGER CSTOP.ONUM,OLAST 

NAMELIST/CASE/HAO,HPO,AO,EC ,XI C ,*C , 00, FO , D ATE, X J20 , R ES, V 
♦ SDEC ♦ SR A, BETA, IPG, I PUT, I CATE, ISTEP, IOCC, IPHOTO, I MANU, IB, 
*ME.FSTOP,OSTOP,TSTOP.DELT,DELF,RS,U,PHI1,PHI2* JD,0VLAP,C 
♦.FPHQT01, F PHOT 02 »L ATI ,LAT2,PLAT,LAMBDA,HAREF,HAD0T»JDREF 
ANGLE (X) = A PCD (X, 360. ) +1 80 . -S I GN ( 180 . , X ) 

C 

1 CONTINUE 
C 

CALL STCCASE 
C 

READ I 5 .CASE ) 

C 

GO TO (2. 3, A), I PUT 
C v 

2 A= AO $ E = FO $ XI = XIO $ W=WO $ 0«C0 $ F=FO 
GO TC 5 

C 

3 AO=(HPO+HAC+2.*RS)/2. 

E0=1 •— IRS+HPOI/AO 
GO TO 2 

C 

4 CALL ORBIT (SOEC.SRA, VINE, I- AO,HPC,EETA, A,E,XI ,W,0,PDEC ,PRA,U,RS) 

F = 0 . $ F 0=0 . 

5 IFl A.LT.O > GO TO 1 
PERI0D=2.*PI*S0RT( A**3/U) 

ANGRAT=2 • *P I /PERIOD 


INF, 

I OL AP , I F R A 
ANGF,C ANGS 
.ISINGLE 
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PANDG 


C 

6 GO Tn(7,F >. IDATE 
C 

7 CALL CAL JlJl (WJC.FJDtWNDfFOf DATE ) 

JO= W JD + F JO 

CO TO 9 

r. 

8 CONTINUE 

WJD=FL GAT ( I F I X I JD>> 

FJD=JD-Wjn 

CALL JIJLC AL (DATF#WJC, FJCtO ) 

C 

9 JDO-JC 

HA = ANGLF<FARFF+H ADOT*( JOJCREF) ) 

C 

CALL VECTOR ( JD.DLCS.P AS , OECE, RAE f CECC,R AC . S X ♦ S Y . SZ • E X ,E Y *E Z tCX f CY t 
*CZ.4) 

C 

CALL OPRINT(C) 

CALL LOGIC(G»TSTEPtFSTEPfJJ ) 

10 CALL LOGIC( 1 , T STt P , F STE P , J J* 

C 

CALL UPCATP(TSTFP f FSTEP) 

C 

IF< < JC- JDn)*86^0C.GT.TSTCP) GO TC I 
IFtONUM.GT.CSTOP) GO TO 1 
IFIFTOTAL.GT.FSTOPI GO TC 1 
C 

IF(OLAST.NE.OKUM) CALL CPP INT ( 1 ) 

C 

GO TO (11.12»13I*JJ 
C 

11 CALL GPP I NT 
GD TO 10 

C 

12 CALL FPRTKT 
GO TO 10 

C 

13 UR I TF ( 6 .101 > 

101 FORMAT ( *0N0 SOLUTION FOUND *) 

GO TO 1 
FND 
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SUBROUTINE PR AT I Cl FOOTO,ECCTF , RATIO ) 

THIS SUBROUTINE CALCULATES THE CVERLAP RATIO BETWEEN TWO FOOTPRINTS 

FOOTO - MATRIX OF 4 COLUMN VECTCRS DEFINING PREVIOUS FOOTPRINT 
FOOT F - MATRIX OF 4 COLUMN VECTCRS DEFINING PRESENT FOOTPRINT 
RATIO - FOOTPRINT CVERLAP RATIC 

CGMMCN/8LCCK3/Pn*DR*UfRStPI»HAREF»HADCTfJDREF 

COMMON/BLOCK5/IPG.IPUT, I DATE, I STEP, I OCC , I PHOTO , I M ANU , IB, IOLAP, 

1 IFRAMEt! SINGLE 
C 

REAL MAGV t NORM* NCRMO * NC PM F , NF JO I , NFLO I , MAG , M AG J I , MAGL 1 1 NO JF I ,NQLF I 
C 

DIMENSION FOOT (3, 6) , FOOTO (3 ,4) , f CCTF < 3 ,4 ) • NORM ( 3 , 6 1 , NCRMO( 3,4) ,NCF 
♦ MFC 3,4) , VI NT (3) ,ANG0(4) ,NFJ0I (3 ) , NFLC I (3) , ANGF (4 I , NO J F i ( 3 ) » NOL F I ( 3 
♦ ) »C AND (4) • C ANF ( 4 ) 

DIMFNSICN C(3) ,X(3) , «Z(8) , I STORE (8) 

ANGLE IX )=AM0D(X,36C. ) +1 8C SI GN (1 60. , X) 

C 

R AT 10=0 • $ EXCFS=0. 

C 

N = 0 
C 

DO 3 1=1,4 

J=I+1 $ IFIJ.FQ.5) J=1 

NOR MO (1,1 l = FCOTO (2 1 1 ) ♦FOOT C I 3 , J )- FOOTO ( 2 , J ) ♦FCQTO ( 3, I > 

NOR MO I 2, I ) = FO0TO(I,J)*FC0TO(3tn~FOOTC(lfI ) ♦FCOTGt 3 , J ) 

NORMOH, I ) = FCCTO(lt I ) ♦FCOTC ( 2 » J )- FOOT 0 ( I , J ) ♦FOOTO (2,1) 

NORMF (1,1 )=FOOTF(2fI)^FOOTF(3tJ)-FOCTF(2fJ)*FOCTF(3 ,1 ) 

N0RMFI2* I) = F0GTF(1,J) ♦FOOT F ( 3 , I )- FOOT F I 1 , I ) ♦FOOTF ( 3 , J ) 

NDRMFI 3* n=FOOTF(l ,1 )^FCCTF(2 , J)-FOCTF ( 1 , J ) ♦FGGTF (2 , I) 

C 

TEMPO =0. $ TEMPF=0. 

DO 1 K— 1 • 3 

TEMPO = TE MPC + NORMG C K , I ) *+2 
1 TEMPF=TEMPF4N0PMFlK f I )**2 
TEMFC=SCPT (TEMPO) 

TEMPF=SQRT(TEMPF) 

IF I ABS(TEMPC).LT..01)G0 TC 8C0 
IF C ABS(TFMPF).LT..01)GC TC SCO 
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PRATIO 


CO 2 K= 1 * 3 

NOR MC ( K » I )=NCRMG(K,I J/TEMFC 

2 NOPMFIK. I)=NCRMF(K, U/TEMPF 
C 

ANGO f I ) = ACC S ( (FOOTOl 1 , 1 MFCOTQ ( 1 , J ) +FCCTOI 2 , 1 » *FOOTO ( 2 . J MFOQTC ( 3 , 
*IMF00T0<3,J> l/RS/RS ) *RC 

3 ANGF ( I ) =ACCS ( ( FGCTFI 1 , 1 ) *FCOTF ( I , J) + FOOTF( 2 » I> *FOOTF ( 2 . J MFOCTF ( 3 , 
♦ IMFCCTFI3, J) >/RS/RS)*RC 

C 

ANG=0. t SUP=0. 

J=4 

DO 4 1 = 1,4 

rANO(I)=ACOS(-l.*{NDRMG(I»J)*NCFPCIl« I) +NORHO ( 2, J J *NORMOI 2 , I) +NGR ¥ 
*013. J MNORMOO, I J ))*RD 
ANG = ANG+C ANO ( I ) 

C4NF { I MACQSI — 1 . * l NORMF if»JMNCPPFIl»I) +NQR MF ( 2 , «H *NOFMF 12,11+NCPP 
*FC3 t J»*NCRNF(3.im*R0 
SUM=SUN+CANFJ I » 

4 J=I 

C 

EXCESC=4NG-360. 

FXCESF=SUM-36C. 

Ml=-1$M2=1$N3=2 

C 

DO 5 1-1 ,4 

K=I+1 $ IF(K.EQ.5I K=1 
DO 5 J = 1 « 4 

L = J + 1 * IF IL .FQ.5 J L = 1 
C 

00 5 ISIGN=P1,M2.M3 

VINTI 1 MNnRPC(2. IMNCRMFI 3, JI-NCFPFI2, JMNORMC (3 .1) 

VINT (2MNCRMFI 1, JMNORMOI3, I I-NCPMOl l, I MNORMFI 3. J» 

VINTI 3) =NOR MO ( 1 » I) *NORMF ( 2 , J >-NC RNF ( I ♦ J ) *NCRMC 12 . 1 1 
MAGV=SQRT (VINK 1 )**2*VINT<2 > 4 *2 < V INTI 3M*2 MI SIGN 

VI NT ( 1) = V I NT 1 1 ) /M AGV % VI NT I 2 ) * V 1 NT 1 2 ) / M4GV $ VINT (3MVINT I 3> /M4GV 
C 

ANGVI =ACOS l (FCCTCC1, IMVINTI1 MF00TCI2, I » *V I NT ( 2 ) +FOOTO ( 3, I MV I 

*NT ( 3 ) MRS MR D 

ANG VK = 4CCS C I FOOTO ( 1 , K M V I NT ( 1 MFOOTO < 2 ,K > *V INT I 2 MFOOTO ( 3, K M V I 

*NT ( 3 ) ) /RSMRD 
C 

IF! 43 S ( ANGVM4NGVK— 4NGCII )).GT..CC1) GC TO 5 
C 
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PRAT10 


C 

C 

C 


C 


c 


c 


c 


c 

c 


ANGVJ=AC(1S( 
*NT(3> >/RS)*RD 
ANGVL=ACOS( 
*NT( 3 ) )/RS)*RO 


( FOOTF ( 1,J)*VINT(1)+F00TF(2, Jl * VI NT ( 2 ) +FOOTF ( 3 , J > * V I 
( FOOTF (l,D*VINT(l) + FOOTF(2,U*VINT(2> +FOOTF ( 3 , L » * V I 


IF( AfiS( ANGVJ+ANGVL-ANGF( J H.CT.l .E-6 > GC TO 5 

TEMP=VINT(1)*FOOTO(1, I) -»V I NT ( 2 ) * FOOTO ( 2, I > +VINT( 3>*FOOTO(3,I> 
IF( TEMP. LT. O.l 41,42 

41 VINT( 1 I=-VINT(1) * V INT( 2 >=- VINT (21 * VINT( 3)=-VINT( 3) 


42 N=N+ 1 

FOOT ( 1 »N ) * VI NT ( 1 1 
FOOT ( 2 ,N ) =VINT( 2) 
FOOT ( 3 » N I =V I NT <31 


5 CONTINUE 


RO B 1*1,4 
SUM=0. 

DO 6 J = 1 .4 

L = J+l $ IFIL.E0.5I L=1 

NF JO I (1 )= FOOTF ( 2 « J ) *FCOTC ( 3 , I I-FCCTOl 2,1 )*FOOTF( 3,J» 

NFJOI ( ? > = FOOTO( 1 ,I»*FOOTF(3» JI-FCCTF ( 1 , J ) *FCCT0 ( 3 , I) 

NFJOI (3 » = FCCTF (1 , J)*F00TC(2, I)-FCOTO( 1 , I) *FOOTF ( 2 , J I 
MAG J I * SOP T ( NFJOI (1 J**2+NF JOI (2)**2+NFJCI(3)**2) 

NFLOI (l» = FOOTF(2.L)*FOOTOl 3, I )-FCOTO( 2, I)*FOOTF( 3,L» 

NFLOI ( 21 =F00TC(1 , I»*FOOTF (3 ,LI-FCCTF ( 1 , L) *FCCTO ( 3 , I ) 

NFLOI ( 3»=F00TF( 1 ,L)*FOOTC( 2 ,1 )-FCCTO( 1 , II *FOOTF ( 2 ,L ) 

M AGL I = SORT (NFLOI ( 1 1 **2+NFL C I ( 2 ) **2+NFL0 1(3 ) ** 2 I 

6 SUM=SUM+ ACCS ( ( NF JO I (1)*NFLC I (1 UNF JC I ( 2 > *NFLO I ( 2 1 +NF JOI ( 3I+NFLC 

* 1(311 /MAGJI/MAGLI )*RD«-ACOS( ( NF JC I ( 1 1 *NCP MF ( 1 , J ) +NF JCI ( 2 ) *NC 

* RMF(2,J) + NFJ0I(3I +NOPMF ( 3 » J ) J/MAGJI !*RD+4C0S(-1*( NFLOI ( 1I*N0RM 

* F( 1,JH-NFL0I<2)*NCRMF(2 , J > + NFLGI (3) *NOR MF ( 3 , JI > / M AGL I ) *R0-180 . 

I F ( ABS( SUM-EXCESP) .LT..0C1) 7,8 

7 N=N+1 

FOOT (1,N)*FCCT0I1»II 
FOOT (2,N1*F0CTQ(2 ,11 
FOOT! 3,N)*FOOTO( 3,1 ) 
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PRATIO 

THE T A I =C AND ( I) 

C 

8 CONTINUE 
C 

no ii 1 = 1,4 

StJM = 0. 
nn 9 j= l *4 

L = J+ 1 i IFtL.EQ.*!) L = 1 
C 

NOJFI (1)=F00TCI2,J)*F00TF(3 , Il-FCCTF ( 2 , I) *FCCTO( 3 , J ) 

NOJFI (2>=FOOTF(1, I)*FOOTO( 3,JI-FCCTO< 1 , J ) *FOOTF ( 3 , I ) 

NOJFI <3)=FCGTCI1 , J ) * FCOTF ( 2 , I )- F CCT F ( 1 , I )*FOCTO( 2 * J ) 
MAGJI=SORT(NOJFIIll**2+NC JFII 2»**2+NOJFI(3l**2) 

NOLFI II) =FOOTCI2,l)*FCOTF<3 ,I)-FCCTF(2, I1*FCCTC(3,L) 

NOLFI I2) = FnOTF( 1, I ) *FOOTO l 3 ,L )~FCCTC I 1 ,L) *FOOTF I 3 , 1 ) 

NOLFI <3)=Fn0T01 1 ,L )*FCOTF(2, I I-FCCTF ( 1, I)*FOOTO(2,L ) 

MAGL I = SQR T ( NOL F 1 1 1 1 **2+NCLF II 2 ) **2+N0LF 1 13) **2 ) 

C 

9 StJM=SUM+ACnS( I NOJFI 1 1)*NCLFI ( 1 1+NOJFI (2 ) *NCLFI 12 ) +NOJF1 ( 3 I *NOLF 

* I (31 J/MAGJI/MAGLI )*PD+ACOSI I NOJF I (1) *NORMOl 1 , J H-NOJF I ( 2 I *NC 

* R NC I 2 »J) + NCJFI ( 3)*NCPPCI3,J I ) /M AG J I) *R O* ACOS (-1* I NOL F I ( 1 1 *NGP F 

* Oil, J)4N0LFI(2)*NQRP0( 2,J)+NCLFII 3I*N0RMCI3, J)l /M6GL I ) *RD-180 . 
C 

IF! ABS1SUP-F XCESOl.LT. 1 . E— € I 1C ,11 
C 

10 N=N+ 1 

F00T(1,M = FlCTF(1,I) 

FOOT ( 2 *NI = FCOTF (2,11 
FOOT I 3,N)=FOOTF( 3, I ) 

THFT A I =C ANF I I ) 

C 

11 CONTINUE 
C 

IFIN.EC.C) GO TO 800 
2 09 NC = N 

00 201 1=1, NC 

PAG=SGRT(FCCTll , 1 ) **2 +FOCT ( 2 , l ) **2«-F0CT(3, I 1**2) 

00 201 J=1 . 3 

201 FOOT I J.I ) =FCOT( J, 11/ MAG 
C 

DO 20 2 1=1,3 

2 02 Cl I)= I FOOT! 1,1) +FCCT ( I,2))/2. 

C 


38 



APPENDIX B 


PRATIO 

MAG=S0RT I C il )**2 + C<2)**2*C(3)**2 ) 

DO 203 1=1,3 

203 CU>=Cm/MAG 

c 

XY=SQRT(Cm**2+C(2)**2) 

ST=C(2)/XY 
CT= C ( 1 ) / XY 
SP=C(3) 

CP = SQR T ( 1«-SP*SP) 

DO 205 1=1. NC 
DO 204 J = 1 , 3 

204 X ( J ) =FCOT ( J , I ) — C f J 1 
C 

YP = -X(ll*ST+xm*CT 

7P=-X (1 l*SP*CT-X(2)*SP*ST+X <3)*CF 

A 7 ( I ) = ANGLE ( ATAN2(YP.ZP ) *RC ) 

205 I STORE ( 1 1 = I 

c 

DO 206 1=1. NC 
M=NC— I 

CO 206 jl = 1 . M 

IF f AZ( JI.LT. A7( J + D) GO TC 2C6 
AT EMP = flZ ( J ) 

I TEMP = I ST ORE ( J) 

AZf J)=A7l J + U 
ISTGRE(J)=ISTCRE(J+1) 

AZ( J+1I=ATEMP 
I ST ORE ( J + l ) = ITEMP 

206 CONTINUE 
C 

DO 2C7 K=1 . NC 
1= I STORE ( K ) 

J= I STORE ( K + l ) 

IF ( K % EQ« NC ) J= I STORE ( 1 ) 

K0RM(1*KI = FCCTI2,1 1 * FOOT ( 3 » J )-FCCT C2.J>*FOOT(3» I) 

NORMI2.K)=FOCT;( 1, J)*FOOTI 3, n-FCCT(l,I)*E0GT(3t Jl 
K0RMI3.KI*F0Ctll» I)*F00T(2r J )-FCCT < 1 , J > ♦FOOT < 2 ♦ II 
MAG = SORT< NORM( 1 ,K)**2+NCRM( 2»K)**2 + NCRV!I3.K 1**2) 

DO 207 L=1 ♦ 3 

2C7 NORM(L.K)=NORW(L.K)/WAG 
C 

J=NC % SUR=0. 

CO 208 1=1, NC 

SUM= SUM* ACOSf-l.M NORM! 1,1 !*NCRR (1, Ji+NCRM(2,I >*N0RM(2, J)+NCRM<3. I 
♦ >*NGRM ( 3 , J ) ) >*RD 
208 J=T 
C 

EXCES=SUM-(NC-?)*180. 

RATIO=EXCES/EXCESO 

c 

800 CONTINUE 
RETURN 
END 
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SUBROUTINE LCG I C IK , T STEP , FSTEP , J J > 

C 

C THIS SUBROUTINE EXAMINES PROGRAM CPTIONS TO OETERMI NE WHETHER THE 

C NEXT EVENT IS A GRCUNDTR ACK CR FHCTOGR APHY EVENT AND CALCULATES 

C THE TIME INCREMENT AND TRUE ANGMALY INCREMENT TO THE NEXT EVENT 

C 

C K=0 INITIALIZE 

C K=1 COMPUTE NFXT STEP ( TSTEP IN SEC AND ESTEP IN DEG I 

C TSTEP - TIME INCREMENT TO NEXT EVENT 

C FSTEP - TRUE ANOMALY INCREMENT TO NEXT EVENT 

C JJ=1 UPDATE AND CALL SUB GPPINT 

C JJ=2 UPDATE AND CALL Sue PPFINT 

C JJ=3 NO SOLUTION 

C 

r.OMMGN/BLCCKl/ JDC. JD , D ATE , CNUM , M , OL AST , CEL T , DELF 
COMMON/BLOCK?/ A, F.XI , W, 0,F , PER ICC, ANGR AT ,F 0 , FTOT AL » TTCT AL 
C0MM0N/BLCCK3/RD»CP.U,RS» PI »HAREFfHADOT. JOREF 
COMMON/eiCCK4/SX,SY, SZ.EX ,EY,EZ ,CX,CY,CZ 

COMMON/R LOCKS/ I PG . I PUT , I C AT E , ISTEP. IOCC , I PHOTO , 1 M ANU , IB . IOL AP , IFR A 
*ME» ISINGLF 

COMMCN/BLfCK6/TSTOP,FSTCP,CSTQP 

COMMON/BLOCK7/HAC,HPO,AO,EC , XI C ,WC . GO . V I NF , SDEC , SR At BETA 
CnMMr.N/BlCCK8/XJ20,RESt PHIItPHI2 ,0 VL AP ,C ANGF ,C ANGS.F PHOTOl f FPH0T02 
♦ .FOOTL tL AMEC A.LAT1 ,LAT2 tPLAT 
COMMON/BLCCK9/NUMP 
RF AL LAT1.LAT2. LAMBDA, LATLAST 
REAL M2.J0TEMP 
REAL JD.JCCtJDREF 
INTEGER CSTCPtONUM.OLAST 
INTEGER FL AGP 

DIMENSION CATEI6) .FOOTLI 3, A) 

CIMENSICN R PCW ( 3 , 3 ) 

ANGLE (X) =AMODIX,3fcO. H-18C.-SIGM1E0. , X) 

C 

C 

KK = K +1 

GD TO II .101) , KK 
C 

1 FL AGP = D 
TNFXTP=0. 

TNEXT G = 0 . 

NUMPsD 
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LOGIC 


GO TO 8CC 


101 TF I ABS(FTCTAL) .GT..001) GC TO 2 
GO TO ( 1C7 * 102 , 1 C2 ) » I PG 

102 GO TO (1C4, 103,105,108) ,IFFAME 

103 CALL RAY(C. ,0.,XC,YC,ZC,KK) 

IF1KK.EQ.C! GG TO 2 

CALL DCT( SX.SY.SZ,XC ,YC,ZC,G) 

IF ( G. LT. PHI2 * AND. G.GT. PH II ) 104,106 

104 FL AGP = 1 
MJM P= 1 
CALL PPRINT 
GO TO 2 

105 DFLF12=ANGLE(FPHCTC2-FPHCTC1 ) 

CEL F02= ANGLE (FPHDT02-F) 

IF( DFLFO? .LT .DELF12 ) 104,106 

106 GO Ttf (2*2,1071 ,IPG 

108 CALL RAY(C. ,0. ,XC,YC ,ZC ,KK) 

IF(KK.EQ.O) GG TO 2 

SI NC = ZC /R S 

IF(LAT1 .GT.LAT2 ) 109,110 

109 TFMP = L AT2 
L AT 2=L AT 1 
L ATI =TEMP 

110 SIN1=SIN(LAT1*DR> 

S I N2 = S I N (1 AT2*DR) 

KK= 1 

IF( SINC.L T. SIN2.AND.SINC.GT.SIN1) KK = 2 
GO TO (106, 1041, KK 

107 CALL GPRINT 

2 GO TO (3,17,8) , IPG 

3 JJ=1 

GO TO (4 ,6) , I STFP 

4 T STFP=DE! T 

5 CALL TCGNIC(U,F, A,F,T1) 

R2= ANGRATMTHTSTEPI 
CALL TINVS(M2,E,EC2,F2) 

F2=ANGLE(F2*R0I 

F STE P-F2-F 

IF ( FSTEP .LT *0. ) FSTEP = 36C . +FSTEP 
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LOGIC 


GO TO ROC 
C 

6 FST FP=DEL f 

7 CALL TCr.UC(U,E.A.F,Tll 

CALL TCCNICUJ.E. A,F*FSTEP ,T2) 

TSTFP = T2— T1 

IF ( TSTEP.LT. 0. ) TSTEP=PEPICOTSTEF 

GO TO ( RPO * 2* * ? £ > ) * I P G 
C 

R CONTINUE 
C 

IFI ARS(TNEXTP) .LT.l.F-b) 17,13 
9 TNEX TP =T STFP 
C 

10 IF(ABSITNEXTGl.LT.l.e-fe) 11,14 
C 

11 GO TO (12.131. ISTEP 
C 

12 TNEXTG = OELT 
GO TO 14 

C 

13 CALL TCCNTCCU.E. A.F.T1 1 
CALL TCCNIUUit, A.F+DELF,T2) 

TNF XTG-T2-T1 

IF( TNEXTG.LT.O. 1 TNE XTG=PE R I GC+ TKE XTG 
C 

14 IFI TNFXTG.LT. TNEXTP1 15,16 
C 

15 TSTE P = TNF XTG 
TNFXTP=TNEXTP-TSTEP 
TNr XTG = 3 • 

JJ=1 

GO TG F 
C 

16 TSTE P« TNF XTP 
TNEXTG=TNEXTG-TSTEP 
TNFX T P=0 • 

FL AGP = 1 
MJMP = NUMP + 1 
J J = 2 
GO TO *5 
C 

17 CONTINUE 
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IF( FLAGP.EQ.O) 29,18 
C 

18 GO TC (18 ,2C. 20,20 .IFRARE 

19 FLAGP=0 
MJMP=0 
GO TO 17 

C 

20 IF( IOLAP.EQ.O) 21,23 
C 

21 GO TO ( 3C7 , 6 I ♦ I STEP 
C 

23 CALL OVERLAP(TSTEP) 

GO TO 308 

307 TSTEP=DEL T 

308 CALL TCCMCdJ, E.A, F,T1 ) 

R2=ANGRAT*( Tl+TSTEP) 

CALL TI\VS(P2.E,EC2,F2) 

F2=ANGLE(F2*RD) 

FSTEP=F2— F 

IFIFSTE’P.LT.C. ) F STEP=36C . ♦ FST EF 
C 

24 GO TO (242,242 ,241,293) .IFRARE 
C 

241 CALL TCCNIC(U.E.A,F,T1) 

K2=ANGRAT*(T1+TSTEP) 

CALL TINVS(M2» E,FC2, FNEXT ) 

FNEXT=ANGLE(FNEXT*RD) 

IF( FNEXT .GT .F PHOT 02) 261,25 
C 

242 JOTFPP= JD+TSTFP/6640G, 

CALL VECTOR ( JOTERP, 81 ,B2,P3 ,E4,e' »B6»$XTEMP»SYTEMP»SZTEMP, 87, 69,68 
*,810. 811,812,4) 

F =F +F STE P 

CALL RAY(0.,0.,XC,YC,ZC,KK) 

F=F-FSTEP 

CALL OCT( 5XTEPP.SYTEyp,SZTEHP,XC ,YC, ZC.GTEMP) 

C 

IF(GTFMP.LT.PHI?.AND.GTEMP.GT. FH 1 1 ) 25,27 
C 

25 GO TO (800.26,9) ,IPG 

26 J J=2 
MJMP=NUMP+1 
GO TO 5 
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27 CONTINUE 
00 28 1=1.2 
DO 23 J= 1 . 2 
F=F+FSTFP 

CALL RAY ( (-1. 1 ** I *C ANGF t (—1# ) * * J *C ANG $ » XC . YC . ZC . KK J 
F=F-F STF P 

IF(KK.EO.O) GO TO 28 

CALL OCT (SXTEKP.SYTEMP, SZTE*F,XC ,YC. ZC.GTEMP ) 
iF<GTEMP.LT.PHI?.AND.GTEMP.GT.PHm GO TO 25 
23 CONTI NUF 

231 FL AG P = 0 
NIJMP =Q 
GO TO 17 

293 F=F+F STEP 

CALL RAY(C. .C. . X T . YT ♦ Z T ,KK ) 

F=F-FSTEP 

IFtKK.EQ.QI GC TC 799 

SINC=ZT/RS 

PLATC = ASIMSINCI*RD 

IF( S I NC • GT. S IN2 • OR. SINC.LT. 5IMI GO TC 296 

294 GO TO 25 

296 F=F+FSTFP 
I ABOVE = 0 
IBFL 1W=0 
DO 299 1=1.2 
no 299 J = 1 • 2 

CAL L RAY( (-1 . )**I*CANGF f (-1. ) ** J *C ANGS , XT . YT.ZT.KK) 
IFCKK.EO.C) GO TO 299 
SLAT=ZT / R S 
PL AT 1= A S I N ( SL AT 1 *R0 
IF! SL AT.GT.SIN2 ) I A BO VE= I A B CVE* 1 
IF( SLAT.LT.STN1> I BE LOW* I BELOW + 1 
799 CCNTINIJE 

IF ( I ABOVE. EG). 4. OR. I BELOW. EC. 4) 2CC.3C5 
300 CENT I NUF 
F = F-F STEP 
GO TO 281 
305 F=F-F STEP 
GO TO 25 
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C 

29 CONTINUE 
C 

60 TO <38,30, 36.521 ,IFRAP£ 

C 

30 CALL SUNBANC(SX.SY,SZfPHHtXI,W»C»RPCW*Fl,F2» I T YPE1 » I7YPE2, 
CALL SUNR AND< SX , SY, SZ ,PHI 2 ,XI , M ,0 ,RPQW,F3 , F4 , ITYPE3 . I TYPE4 , 
IF< ITYPE1.EQ.0) F1 = IC'0C • 

IF< TTYPE2.E0.01 F2=1000. 

IF I ITYPF3.EC.0 1 F3=1000. 

IF( I TYPE4.EQ.0) F4=100C « 

IF(AMIN1<F1,F2,F3,F4).GT. 999.1 GC TO 799 

0F=AMIN1< ANGLE (Fl-F). ANGLE < F2-F 1 , «NGLE < F3-F 1 , ANGLE (F4-F 1 1 
FSTEP=OF 

IF C ABS(r)F-ANGLFfFl-F) 1.LT..CCI1 GC TO 31 
IF < ABS ( OF— ANGL E < F 2— F 11.LT..CC1) GC TO 31 
GNEXT = PHI 2 
GO TO* 32 

31 GNEXT = PH 1 1 

32 CONTINUE 
C 

IF ( I PHOTO . E 0,0 1 GO TC 37 
C 

NC N T = 0 

CALL TCONICIU.E.A.F.Tl 1 
CALL TC0NIC(U.F.A,F+DF,T21 
DT=T2-T1 

IF(0T.LT.0.1 DT=PER I CO+CT 
C 

CALL UPC ATE ( CT , DF 1 
T1 = T2 

CALL RAY (0.,0.,XC,YC, ZC,KK I 
CALL OCTI £X,’SY.SZ,XC,YC,ZC,G11 
C 

33 CONTINUE 

CALL TC0NIC(U.E,A,F*1.,T21 
TSTFP=T2-Tl 

IFITSTEP.LT.O.) TSTEP=PER 1 00 +T STEP 
nF=OF+l. 

OT=DT +TSTEP 

CALL UPDATEfTSTEP.l. 1 

CALI R AY < 0 • , 0 • , XC , YC . ZC » KK ) 

CALL CCTtSX ,SY.SZ,XC,YC,ZC,G21 
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PGwRTF= ( G2— G1 ) / 1 . 

FSTEP=-(G1-GNEXT) /PGWRTF-i. 

T1 = T2 

CALL TCCNIC(U.Ff A.F+FSTEP ,T2) 

TSTEP=T2-T1 

I F f F STF P • GT • O* . AND. T STEP. LT .0 . ) TSTEP=PER IOD+TSTEP 
IF! FST^P.LT.G. .AND.TSTEP.GT .0. ) T $TEP*T STEP-PERI 00 
0F= DF + F ST E P 
DT=DT +TSTFP 

CALL l)PCATE(TSTEPt FSTEP ) 

CALL RAY(0. ,C. tXCtYC .ZC ,KK) 

CALL DOT< SXtSY,SZ»XC,VC,ZC,GI> 

Tl = T2 

IF ( Afi S ( Gl-GNEXT I .LT..001) 25,34 

34 NCNT=NCNT + 1 

IF ( NC NT.FQ.10) 799,33 

35 CONTINUE 
T STEP = D T 
FST F P = DF 

CALL UPDATEI-OTt-DF) 

GO TO 50 

36 CONTINUE 

FSTEP^ANGLE(FPHOTOl-F) 

37 CALt TCCNICf UtEt A.F.Tli 

CALL TCGNIC(U,E*A,F+FSTEP,T2) 

T STE P= T2 -T1 

IF( TSTEP.LT.O . ) TSTFP=TSTEP+PERICC 

GO TO 50 

3B GO TO ( 39 *46 I ♦ T S INGLE 
39 SCI = S INI PLAT*DR) /SIN(XI*CR) 

IFC ABS(SCl) .GE.l .1 GO TO 7^9 
Cl = A S I N ( SCI ) *R D 
C2 = S I GN I 180 • • Cl ) - Cl 
C0 = ANGL E ( W+F ) 

IFC ANGLE (Cl-CO).LT.. 001 ICl=CG-.l 
IF( ANGLE (C2-C0) .LT.. 001 )C2=CC-.l 
CELOAMIN1I ANGL FI Cl-CO) , ANGLE IC2-C0 ) ) 

41 CALL TCCNICIUtEt A.F.TII 
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CALL TCCNICIU.E, A,F+CELC,T2 > 

DT=T2-T1 

IFtOT.LT.O.) DT=PER IOD+CT 
OST AR=0 
WST AR = W 

CALL hELlSC?S,U.XJ20, A, E , X I ST A F , OST «R, DT ) 
FL = ANGLE ( CC+DELC— WST AR > 

FSTEP=ANGLE(FL-F) 

TSTEP=DT 

IF € IPHOTO.EO.OI GO TO 50 

42 NCNT =0 

CALL UPCATEITSTEP.FSTEP) 

CALL RAY(0.»0.*XT»YT • ZT »KK ) 

IF(KK.EO.C) GO TO 799 
PL AT1 = AS1N(ZT/RS )*RD 

43 T1=T2 

CALL TCON IC (U»E»A»F+1 »»T2 I 
0T=T2-T1 

IFJOT.LT.O.) DT=PER I OD+DT 
CALI UPDATE COT , 1 • ) 

TSTEP=TSTEP*DT $ FSTEP=FSTEP+1. 

CALL RAY(0..0.,XT,YT,ZT,KK) 

IF(KK.EQ.O» GO TO 799 
PLAT2=ASIN( ZT/RS »*RD 
PLWRTF=(PLAT2— PLAT1I/1. 

0F= C PL AT— PL ATI I / PLWRTF— 1 • 

T1 = T2 

CALL TCCNICLU.E. A,F+DF,T2 > 

DT=T2— T1 

IFIDT.LT.C. .AND.OF.GT.G .) CT=PEPICO+DT 
IFIOT.GT.O. .AND.DF.LT.O.) CT=CT-FERIGO 
CALL UPC ATE ( CT . OF ) 

TSTEP=TSTEP+DT « FSTEP=FSTEP+DF 
CALL RAYCO.,O..XT,YT, ZT,KK ) 

IFLKK.EO.Ol GO TO 799 
PL AT1 = ASIN( ZT/RS )*RD 
IF(ABS(PLAT1— PLATI.LT..CC1 ) 45,4* 

44 NCNT=NCNT +1 

IF( NCNT .F0.10 I 799,43 

45 CALL UPCATcL— TSTEP,— FSTEP) 

GO TC 50 

46 CONTINUE 
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DUMMY LOGIC 


ESTEP=10C. 

GO TO 37 
C 

52 SXI=SIN(XI*DR) 

SCI-$IN(LAT1*DR) /SXI 
SC3*SIN(L AT2*DR)/SXI 
C0 = ANG L E ( W+F ) 

I F ( AB S ( SC 1 ) *GF • 1 • ) 55,56 

55 C1=C0+359.C9 % C2=CI $ GC TC 57 

56 CONTINUE 
C1=ASIN( SCI 1*PD 
C2-S I GNC 1 BC • *C 1 ) — C l 

57 IF( AB S t SC3 ) «GE • 1 • ) 58,59 

58 C3=CO+359 *99 $ C4=C3 $ GC TO 6C 

59 CONTINUE 

C3= A S I N ( SC3 ) *R D 
C4=S I GN ( 1 8C • § C3 )— C3 

60 CONTINUE 

CMIN=AMIN1(C1,C2,C3,C4) 

IFCCMIN.GT.099 , ) GO TO 799 

DEL C= AM INI ( ANGLE (Cl-CO) , ANGLE (C2-C0 1 , ANGLE ( C3-C0I , ANGLE (C4-C0) I 

IF( ARSCDELC— ANGLECC1— CC)) .LT..001) GO TO 53 

IF f ABSC DEL C— ANGLE (C2— CO) ) .LT..CC1) GC TO 53 

PL AT = L AT2 

LATLAST=L ATI 

GO TO 54 

53 PL AT = L ATI 
LATL A S T = l AT2 

54 GO TC 41 
C 

50 GO TC (51,51 ,9) , IPG 
C 

51 FI. AGP = 1 
NUMP a 1 

J J=2 

GO TO 800 
C 

799 J J=3 

800 RETURN 
END 
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subroutine gprint 

c 

C THIS SUBROUTINE CALCULATES AND WRITES GROUNDTRACK DATA 

C 

DIMENSION DATE<6) ,OTME<4) ,F0CTL(2 ,4) 

REAL JCC, JD.LAT3.LON2 

REAL LAT1 .LAT2. LAMBDA, LATLAST 

INTEGER CSTCP.CNUM.OLAST 

C0MM0N/RL0CK1/ JDO. JO, DATE . C.NUM.Hf- ,0L AST, CELT »DELF 
COMMCN/BLOCK2/ A, E.XI , W ,0, F , P ER ICC, ANGR AT , FO ,FTOT AL , TTOT AL 
COM MON /BLOCK 3 /RD ,DR , U ,R S , P I ,H AREF ,H ADOT , JDREF 
COMMON/ BLCCKA/SX.SY.SZ, EX, EY, EZ ,CX ,CY,CZ 

COMMCN/BLCCKfi/IPG.IPUT , ICATE, ISTEP, ICCC. I PHOTO, IMANU, IB, IOLAP, 1 FR A 
♦ME. ISINGl F 

COMMO N/B LCCKfc /T STOP,FSTCP,CSTCP 

COMMON/ BL0CK7/HA0.HP0, AO, EC , X 10 , WC , CO . VI NF , SDEC , SR A, BET A 
C0MM0N/BLCCKB/XJ20,RES,PHI1,PH!2,CVLAP,CANGF,CANGS , F PHOTO 1 , FP HOTO 2 
♦, EOOTL. L AMBD A. L ATI, LAT2.PL AT 
ANGLE (XI =AMCD(X,3fcO. I+1BC. -SIGN lieo. , XI 
CALL CCNCAR( A,F ,XI . W ,0 , F , X 1 , Y 1 , Z1 ,DX1 ,0 Y1 , D Z 1 ,U) 

X2 = Xl*C0S(HA*DR »+Yl*S IN ( H A ♦DR I 
Y?=— X1*SIN(HA*DR I+Y1*C0S(HA^DRI 
Z?=71 

0X2=0X1*CCS(HA*DR1+DY1*SIN(HA*CR1 
DY2=— DXI*SIN(HA*DR)+DY1*CGS(KA*DRI 
D72 =D Z1 
C 

CALL DOT (XI, Y1 ,Z1,SX,SY,SZ,PHII 

CALL CARCON(X2.Y2,Z2,DX2,DY2,OZ2,A2,E2 ,X 1 2 , W2 ,C2 , F2 . U I 
C 

CALL CAR SPM ( X2 .Y2.Z2 ,DX2,DY2 ,CZ2 ,R2 ,L0N2 ,LAT3, V2 » GAM2 ,SI G2 I 
L0N2= ANGLE (L0N21 
H2=R2-RS 
C 

WM=HAD0T*DR/B640G. 

VHI=V2*C0S(GAM2*CR) 

VHIE=VHI*SIN(SIG2*0R) 

VHI N= VHI *C0 S( S I G?*DR I 
VMF=RS*COS ( L AT3*DR 1 *WM 
VTE=VHI E-VME 

SIG2R=ATAN2(VTE,VHIN)*R0 
VHR=SQRT( VHIN**2*VTE**2) 
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HVOM = \HR /H2 

CEL jn = TT0T&L/86*0C. 

DAYS=DEl JC- A MOO ( OF L J D »1 . ) 

DELHRS=<CFL JD-DAYS)*24. 

HRS=DELHRS-AMOO(DELHRS f 1. ) 

CEL M IN= ( OFL HRS— HR S ) *60. 

KIN=DFLM! N- AMOf) (DELMIN»1.) 

SEC=( DELMIN-MIN) *60. 

CT M E ( 1 ) = D A Y S 
OTMFI 2 ) =MR S 
CT M E ( 3 ) = M I N 
0TME(41=SEC 
S=TTOT At 
C 
f. 

WRITF (6.1 Cl) J0,CATE,C1>E,S,CMJR 

WRITE (6. 102)LAT3.10N2»H2.V2,GAM2,PHI ,H VOH .H A , W .0 .F 

101 FORM AT ( 1 HO /* JULIAN D AT E * » F 18 . R , 2 X , *C AL ENO AR CATE*.F4,0.4F3.0 »F1 . 2 
1,2X.*0RBIT TIME*,3F3.0»F6.2»2X,*CSEC1*»G12.6.*0RBIT NC*.I2I 

102 FORMAT (1HC ,*L AT = * » F It- .8 » 2X » *L CNG =* ♦ E16. 8 , 2X , * ALT =*.E16.8»2X, 

1 *VE L I =*.E16.P.2X,*FPAI = * . E 16 . £ , /I X.*PH I = * .E16 .8 . 2 X . *VOHR =*t 
2E16.3,?X,*HA =*.E16.8,2X f *W =* . E 16 .8 . 2 X , *0 = * .E 16. 8 » /I X ,* T 

3. A. =*.816.8) 

RFTIJRN 

END 
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SIJ8R OUT INF CPRINT(K) 

c 

C THIS SLRFOUTINE KRITFS THE PROGRAM INPUTS 

C 

C K=0 I NIT AL PRINTOUT ( START OF CASE I 

C K=1 PFRIAPSIS PRINTOUT 

f. 

CaMM!CN/BLCCKl/JCC,jn,CATE,CNUM,hA.OLAST,OELT,DELF 
COMMnN/RLOCK2/A»E,XI*W,O,F»PERI0D*ANGRAT»FO,FTCTAL*TTOTAL 
COMMON/BL CCK3 /RD *DR*U*RS»PI»HAREF*HADGT*JOREF 
COMMON/ RL CCKi / SX, SY, $ Z, E X, E Y, E Z, C X ,CY,CZ 

COYMON/BICCK5/IPG,IPUT,ICATF, ISTEP, IOCC, IPHOTO* IMANU, IB , IOL AP , I F R A 
*ME. I SINGLE 

COM MC N/ BL CCK6/TSTOP, F STOP, CSTCP 

r.nMMON/RLOCK7/HAC,HPO,AG,En,XIC,fcC ,00. VINF ,SOFC , SR A, BETA 
CCMM0N/RLCCK8/XJ20*RES,PhIl,PHI2,0VLAP*CANGF ,C ANGS,F PHOTO 1 .FPH0TG2 
*.F0nTL,LAMBD&,LATl.lAT2,PLAT 
Cl MENS I CN CATEL6 J.FOCTL (3 ,6 ) 

REAL LAT1 ,LAT2, LAMBDA ,LATLAST 
INTFGFR CSTCP 
C 

KK=K+1 

GO TO (1,2) « KK 
1 CONTINUE 

WRITE(6»lG0i)A0«RETA*CANGF,CANGS*CATF(l)»CATE(2)»DATE(3)*DATE(4)»C 
1 ATE(5)»CATE(6)*DELF,DELT,EC,FC,FPFOTni»FPHOTC2,FSTOP»HAG,HPO,IB,IC 
2ATF , I FRA ME , I MANU , I OCC , I CL AF , I PG , I PHOTO , I PUT , I S I NGL E, ISTEP, JD, 

3L AMBD A * L ATI ,L AT2 ,00 , CSTGP , GVL A P , PR 1 1 , PHI 2, PL AT , RE S, R S , SOEC , SR A , 
4TSTGP»U»VINF , V«p , X 1 0 , X J20 

1001 FORMAT (*1 LEASE*//* AO =*,Elf.8,* BETA =*,E16.8,* 

1CANGF =*,E16.8,* CANGS =*,616.8/* C ATE C 1 ) =*,E16.8,* 

20 ATE ( 2 ) =*,F16.8,* 0ATE(3) =*»E16.8,* 0ATfc(4) =*,E16.9/* 

3 DATE(S) =*,F16.P.* CATE(o) =*,E16.8,* DELF =*,616.8,* 

4 DEL T =* ,E16. 3/* EC =*.E16.8,* FC =*,E16.3,* 

5 FPHCT01 =*»F16.3» * FPH0T02 =*,EI6.3/* FSTOP =*,E16.6, 

6* (AO =*»E16.8»* HFC =*»E16.8,* IB = *,I 

716/* ICATE =*,116,* I FRAME =*.116,* IMANU =*,116,* 

8 IOCC =*,116/* IOLAP =*,116,* IPG =*,116,* IF 

9H0T0 =*,116,* I PUT =*,116/* ISINGLE =*,116,* ISTEP 

*=*,116,* * , 16X , * JC =*,116/* LAMBDA =*,E1 

*6. B * * L A T 1 =*,E16.e,* LAT2 =*,E16.8,* 00 

*=* » FI 6 . 3/* CSTOP =*,116,* CVLAP =*»E16.8,* PHI1 =* 
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OPRINT 


*.E16. 

8,* 

PH I 2 

= *,Elfc.8/* PL AT 

=*tE16.8t* RES 

= * 

*.F16. 

8.* 

RS 

= * t E 1 6 . 9 t * 

SDEC =*.E16.8/* SR A 

e 


.3.* 

T STOP 

=**E16.8»* 

U = * » El 6 • 8 ♦* 

VINF 

* = 

*♦816. S/* WO 

=* ttl6.fi t* 

XIO = *• E16 . 8# * 

XJ2C 


* = * • E 1 6 • 8 I 

IFdPG.EO.il WRITF<6»106> 

106 FORMAT (*CGP0l)ND TRACK OPTION*) 

IF f IPG.EC.21WRITEC6.107I 

107 FORMAT (*OPHCTO OPTION * ) 

IF("lPG.E0.3IWPITE.(6 #108) 

103 FORMAT ( *0 GROUND TRACK ANC PHCTO CPTICN*) 

IF( IPLT.FQ.l ) WRITE <6. 10 9) AC, EC 

109 FORMAT!* INPUT *,* AO =*,E16.g f * EC =*,E16.8) 

IF! I PUT. EC. 2 ) WRITE (6*110 )FPC, HAG 

110 FORMAT!* INPUT * ♦ * HPC =*,E16.E.* HAC =*,E16.8) 
IFdPUT.F0.31WRITFI6.nl) 

111 FORMAT!* ? VECTOR INPUT*) 

IF! ICATE.EO.l ) WRITE (6 , 1 12 1 

112 FORMAT!* CALENDAR DATE INPUT*) 

IF! I CATE. EC. 2) WRITE (6,112 ) 

113 FORMAT!* JULIAN CATE INPUT*) 

IFdSTEP.FC.l ) WRITE (6,1 14) CELT 

114 FORMAT!* TIME STEP OF *,E16.8,* EMPLOYED*) 

IF! ISTEP.F0.2 )WRITE(6, 115)CELF 

115 FORMAT!* TRUE ANCMALY STEP CF *.E16.8.* EMPLOYED*) 
IF! IOCC.EO.O ) WRITE! 6.1161 

116 FORMAT!* NC CCCULATICN CATA*) 
IF(I0CC.FQ.1)WRITE(6*117) 

117 FORMAT!* CCCULATICN DATA GIVEN*) 

IF! IPHOTG.FO.O) WRITE! 6.1181 

118 FORMAT!* VFRTICAL PHOTOGRAPHY*) 

IF! IPH0TC.E0.1)WRITE(6,119) 

119 FORMAT!* NCN VERTICAL PHOTOGRAPHY*) 
IFITMAMJ.EG.OWRITE (6.120) 

120 FORMAT!* NO MANEUVERS MACE*) 

IFdMANU.EO.l ) WRITEC6.121) 

121 FORMAT!* MANEUVERS COMPUTED*) 

IFdR.EQ.ClWRITE (6.122) 

122 FORMAT!* RMATRIX NOT COMPUTED*) 

IF! IR.EQ.l >WRITE <6 .123) 

123 FORMAT!* RMATRIX COMPUTED*) 

IF!ICLAP.EC.0)WRITE(6 f 124) 

124 FORMAT!* PICTURES TAKEN CN STEP*) 
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OPRINT 


IF! I0LAP.EC.1)WRITE(6 ,125)CVL«P 

125 FORM AT ( * PICTURES TAKEN CN *,E16.8,* OVERLAP*! 

IF! IFRAME.EO.l ) WR ITE ( 6 » 126) 

126 FORMAT!* SINGLE PICTURE TAKEN*) 

IF! IFRAME.EQ.2)WPITE(6,127 ) PH II, PH 2 

127 FORMAT!* PICTURES TAKEN BETWEEN PHI1 OF *.E16.0,* AND PHI2 

10F *,E16.8) 

IF! IFRAME.E0.3)WRITE!6i 128) FPHCTC1 ,FPH0T02 

128 FORMAT!* PICTURES TAKEN BETWEEN FPHOTOl OF *,F7.2,* AND FPHC 
1 T02 OF *,F7.2I 

IF! IFRAME.EO.4J WRITE <6,129 )LAT1,LAT2 

129 FORMAT!* PICTURES TAKEN BETWEEN L ATI OF *,F7.2,* AND LAT2 

1 OF * * F7 . 2 ) 

T0STCP=360.-FG+360.*(0ST0P-1) 

TFSTOP=FSTOP 
B=AMOO!TSTOP, PERIOD) 

CALL TCONIC«U,E, A.FO.TFC) 

XM= ( TFO+B ) * ANGP AT 
IF! XM .GT .2 .*PI ) 3, A 

3 B=B— PER I OD 
XM=XM-2.*PI 

4 CONTINUE 

CALL TINVS!XM,E,EC,FB) 

FOB=FB*RD-FC 

IF!FOB.LT. .001 ) F0B=360 . +FOB 
TTST0P=36C.*( TSTCP-B)/PERICC+FOE 

IF(TOSTOP.LT.TFSTOP.AND.TCSTCP.LT.TTSTOP)WRITEI6,20CO)OSTOP 

2000 FORMAT!* PROGRAM WILL STCP CN OSTOP = *,I5) 

IF!TFSTOP.LT.TTSTOP.AND.TFSTOP.LT.TCSTOP)WRITE!6,2001)FSTOP 

2001 FORMAT!* PROGRAM WILL STCP CN FSTOP = *,E16.8,* DEGREES*) 

IF ( TT STOP. LT. TF STOP. AND.TTS TOP. LT.TOS TOP) WRITE! 6, 2002 )TSTOP 

2002 FORMAT!* PRCGR'AM WILL STOP CN TSTOP = *,E16.B,* SEC*) 

C 

2 CONTINUE 
C 

RETURN 

END 
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SUBROUTINE 0 V E R L AP(TSTEP) 

THIS SUBROUTINE CALCULATES TEE TIRE INCREMENT TC THE NEXT CAMERA 
FOOTPRINT FOR A GIVEN OVERLAP PATIO 

TSTEP - TIME TO NEXT CAMERA FOOTPRINT 


C 


C 

c 


c 

c 


CO^MCN/RLOCKl/JDO, JD, DATE.CNUM.M ,OL AST.DELT ,DELF 
COMMON /BLOCK?/ A, E , XI , W f 0 , F , PERI C C , ANGR AT , FC , FTCT AL. TTOT AL 
CGMMON/fiLOCK3/RO,QR,U,RS,PI,HAREF,HAOOT, JDREF 

C0MM0N/BL0CK3/I PG ♦ I PUT , I C AT E , I STEP, ICCC, I PHOTO, IMANU, IB, IOLAP, IFRA 
♦ ME, IS INGL F 


COMMON/BLCCKB/XJ2C , RES, PH II ,PHI2 ,GVLAP,CANGF,CANGS, 
♦.FOOTL.L am bo A, L ATI.L AT2,PLAT 
REAL L ATI ,L4T2 .LAMBDA, LATLAST 

DIMENSION FOOT 0(3,4 > , FOOT F ( 3,4 ) , K (4 ) *KF (4 ) , C ATE( 6 ) , 

ANGLE (X) = AM0D(X,3SC . I +180. -SIC-NI 180., X) 

CALL RAY (-CANGF, -CANOS, F00T0(1,1 I.FOOTOI 2, 1 I ,F00T0( 
CALL RAYI-CANGF, CANGS»FCCTC(I»2), FOOT 0(2*2), FCOTOI 
CALL RAY! CANGF, C ANGS , FOOTO ( 1, 3 ) , FOOTO ( 2 , 3 I ,FOOTO( 
CALI R4Y( CANGF ,— C ANG S,FCOTC(I ,4 ) , FOOT 0(2, 4) ,FOOTO( 

IF(K(l)+K(2)*K(3)+K(4).EC.4) GO TO 7 

0F=O $ T STF P=0 $ CALL TCCM C ( U , E , A , F, T1 ) 
a TSTEP=TSTFP+1C0. 


FPHQT01.FPHOT02 


FOOTL (3,4) 


3.1) .K(1)) 

3 .2 ) , K( 2) ) 

3.3) ,K( 3) ) 

3.4) ,K(4) ) 


T2=T1+TSTEP 


XM= ANGR AT *T2 

CALL TINVS( XM.E.EC2 ,F2 ) 

FSTFP = ANGLE ( ANGLE (F2*RD)-F) 
DF=OF+FSTEP 

CALL UPDATEd CO. ,FSTEP) 

CALL RAY( -CANGF, -C ANGS, X,Y, Z,K(1 ) ) 
CALL RAY ( — C ANGF , C ANGS , X , Y , Z , K ( 2 ) ) 
CALL R AY ( CANGF, C ANG S , X , Y , Z , K ( 2 ) ) 
CALL R AY ( CANGF, -C ANGS , X , Y , Z , K (4 ) ) 
IF(K(I)+K(2)+K(3)4K(4).NE.4) GO TC 8 
CALL IJPCATEl-TSTEP.-OF) 

GO TO BOO 
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OVERLAP 


7 CONTINUE 
C 

$H4 = $ I N ( H A*DR ) * CMA=CCS(HA*CR) 

C 

CO 1 1=1 

XTEMP = FOGTC( 1 1 T I t Y T£MP=FCCTC ( 2 , I) 
FOQTOtl t I J=XTEMP*CH4+YTEMP*ShA 
1 FOniO(? t I >=-XT£MP*SH4+YTe*P*CH4 


CT =0 • $ OL AP 1 = 1 • 

DF=0. i F S TF P = 1 • 

C 

CALL TCCNIC(U*F»Atf*Tn 
C 

2 CONTINUE 
C 

CALL TCCNICCU.E. A . F+FST £ P , T 2 > 

TSTEP=T2— T1 

IF(TSTFP*LT.O.. AND. FSTE P. GT .O. ) TST£P=P£R IOO+T STEP 
IF(TSTEP.GT.0..4ND.FSTEP.L7.G. ) T S TEP = T STE P-PE R I 00 
OF=DF+FSTEP 
DT=DT+TSTFP 

CALL UPCATEfTSTEPtFSTEP ) 

CALL RAY(-CANGFf-CANGS,FCaTFIi ,1) tF00TF(2*l ) ♦ FOOT F ( 3 . 1 1 * KF (I) I 
CALL R AY (-C ANGF t C ANGS , FCOTF ( 1 , 2 ) , FCOTF ( 2 f 2 ) *FOOTF ( 3 , 2 ) ♦ KF { 2 >) 
CALL RAY{ C ANGF * C ANGS, FOCTF ( 1 , 3 ) t FCOTF i 2 , 3 > tFCOTF ( 3 f 3 1 . KF ( 3 > ) 
CALL R AY ( CANGF*-CANGStFOOTF( 1 * 4 ) ♦ FOOT FI2t4ltF00TF(3»4)*KF{4) I 
C 

SH4 = S I N ( H4*DR ) % CH A-COS ( F A *CR ) 

C 

CO 3 1*1,4 

XT£MP=FfOTF(l #1 ) % Y TEMF = F CCT F ( 2 f I ) 

FOOT F < 1 * I I = XTEMP *CH A+YTE^P*SH A 

3 FOOTF ( 2 ♦ I )*-XTF^P*SHA4*YTEFF*CHA 
C 

CALL PRATICCFCOTCtFOCTF t R AT 10 1 
!F(RATIO.GT..G99S9)GO TC 6 
C 

I F ( ABS (RATIO) .GT **000001 ) GO TO 5 
CALL UPCA7F(~TSTEP*-FSTEP) 

DF = DF-f STEP $ DT = DT-TSTEP 
FSTEP=FSTEP/2. 
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GO TO 2 

5 CONTI NUF 

IF(ABS(RATICHOVLAP).LT..C001) GG TO 4 
POWRTF = (R ATIO-OL API ) /FSTEP 
FSTEP-— (PAT IC— G VL AP ) /POWRTF 
IF( PGURTF.GT.C.O) FSTEP=5. 

IF ( AB S f FSTEP ).GT.*>. ) F STEP* S I GN ( 5 • *F STEP ) 
CLAPI= RATIO 
T1 = T 2 
GO TG 2 

6 F STFP =S . 

T1 = T? 

GO TO 2 

4 CONTINUE 
TST FP = OT 

CALL UPDATE (-CT.t~DF ) 

fl 00 RETURN 
END 
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SUBROUTINE PPPINT 

THUS SUBROUTINE CALCULATES AND WRITES PHOTOGRAPHY DATA 

COMMON/ BLCCK1/ JDC, JD, CAT£,CNUM,E A,OL AST , DEL T, DEL F 
COMMON/BLCCK2/A,E,XI ,W,0,F, PERIOD, ANGR AT ,FO ,FTOT AL ,TT CT AL 
CaMMCN/RLCCK3/RD,DR,U,RS,PI,F4REF,HADOT, JDREF 
C0MM0N/BLCCK4/SX,SY,SZ,EX ,EY,EZ,CX,CY,CZ 

COMMON/ BLOCK5/ IPG, I PUT, IC ATE, ISTEP, IOCC , IPHOTO, I MANU, IB, IOLAP , IFR A 
♦ME , I S I NGLE 

COMMCN/BLCCK6/TSTOP,FSTCP,CSTOP 

COMMON/ ALCCK7/H AC, HPO, AC, EC , X I C, WC , CO , V INF , SOEC.SRA, BETA 
COMMON/BLOCK8/XJ20.RES,PHI1 , PHI 2 ,C VL AP ,C ANGF ,C ANGS ,F PHOTOl ,FPHCTC2 
♦.FOOTL.LAMEC A, L ATI, LAT2.PL #T 
C0MM0N/BL0CK<3/NUMP 
REAL LATI.LAT2, LAMBDA, LATLAST 
DIMENSION DATE I 6 ) ,F00TL (3,4) 

DIMENSION FCCT(3,4),KK(8) 

ANGLE (X>=AMOO( X.36C. ) +1 80 SI GM 1 80 . ,XI 
C 

CALL GPRINT 
C 

CALL RAY (—CANGE »— C ANGS , FC C T ( 1 , 1 ) , FOOT ( 2 , 1 ) , FOOT (3»1).KK(1)I 
CALL RAY (-CANGF, C ANGS , FOCT( 1 , 2 ) .FOOT ( 2 , 2 ) . FOOT ( 3, 2 ) ,KK ( 3 ) I 
CALL RAY ( C ANGF , C ANGS , FC CT U , 3 ) , FOOT ( 2, 3 ) , FOOT < 3 , 3 I ,KK ( 5 > ) 

CALL RAY ( CANGF.—C ANGS, FOOT! 1,4) »F00T(2»4) »FC0T(3,4) ,KK(7I I 
C 

IF(KK(1)+KK(3)+KK(5)+KK(7).NE.4I GO TC 300 
C 

SHA = S I N ( H A*CR ) t CHA=COS ( H A *DR ) 

C 

no l i=i,4 

XTEMP=FOOT (1,1) $ YTEMP=FCCT( 2,1) 

FOOT ( 1 . I ) = XTEMP*CH4+YTENP*SH4 
1 FOOT ( 2,1 ) =-X TEMP* SHA+YTEMP *CH A 
C 

CALL PRATIC(FOOTL, FOOT, RATIO) 

C 

no 201 1=1.3 
DO 201 J= 1 ,4 

201 FOOTH I . J)*FCCT( I , J) 

C 


57 



APPENDIX B 


PPRINT 


WRITE( 6*101) R A T I C » NUMP 

101 FORMAT ( * PHOTO OVERLAP = * F7 • 3 ? 1C X * *F00TPR I NT NUMBER = *13/12X,*LA 
*TI TUBE *.l 1 X, *LONGI TUDE* t 12X ,* SUN ANGLE* t 10 X f *SL ANT RANGE*t7X f *ST AT 
*IC R E SOLUT I CN* ) 

C 

CAl L CCNCARl A*F*XI *W,0,F t X,Y, Z f O>.DY.DZ,U) 

C 

XTEMP = X * YTFMP= Y 
X=XTEMP*CH A+YTEMP*SH6 
Y=- XTFMP^SH A + YTE M P*CH A 
C 

X SIJN -SX*CHA+SY*$HA 
YSUN=-SX*$HA+SY*CHA 
Z SUN = S Z 

c 

N=n \ 

9 N=N+l 

GO TO (ir*lltl2fllfl3»lltl^ 1 1 1 » 15 *200 ) * N 

c 

10 CALL RAY(G..C.fXT*YT,ZT f K> 

GO TO 16 

C 

11 J=N/2 

XT= FOG T 1 1 *J> 

YT = FOO T ( 2 ♦ J ) 

ZT = F OCT ( 3 * J ) 

GO TO 17 
C 

12 CALL RAY(-CANGF*C. *XT f YT ,ZT ,KK< h-1) ) 

GO TO 16 

13 CALL RAYfO. *C ANG S t XT , YT , Z T ,KK ( K- 1 M 

GO TO 16 

14 CALL RAYICANGF.C. . XT * YT t Z T , KM N- 1 )) 

GO TO 16 

15 CALL RAY(0.*-CANGStXT »YT tZTfKKH-1) ) 

GO TO 16 

C 

16 XTF MP = XT * V TF MP = YT 

xt=xtemp*cha+ytemp*sfa 

YT=-XTEMP*SHA*YTEMP*CHA 

C 

17 CALL LATL NG( X T . YT * Z T * XL A 7 ,XLCNC ) 

XLONG = ANGL E ( XL ONG ) 
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PPRINT 


C. 


C 


c 


c 


c 


CALL nOT( XT.YT.ZT.XSUN, Y$ IN , Z SUN ,PH I ) 

$R=SQRT ( ( X-XT)**2MY-YT )**2+(Z-ZT 1**21 
SRE S= SR *R E S 

GO TO <18.19.19. 19,19,19, 19 ,19,19) f N 

18 WR I TE ( 6 * 1 €2 ) XL AT . XL ONG , PH I . SR , S P ES 

102 FORMAT J* CENTERS * 6 X . F6 . 2 * 1 2 X , F7 .2 . 1 3 X ,F7 • 2 * 1 3X . F9 . 3 • 1 3 X , F9 • 2 I 
GO TO 9 

19 I =N— 1 

WR I TE ( 6 . 1 C3 I I .XLAT.XLCNG.PHI ,$R,SRES 

103 FORMAT ( * PC I NT*. 12. 5X.F6 • 2 , 13X, F7 . 2 * 13X.F7. 2 1 13X .F9. 3 *13 X , F9. 2 ) 
GO TO 9 

300 WRITE<6.104) 

104 FORMAT ( *0 PHOTO FOOTPRINT CFF PLANET *) 

200 CONTINUE 


RETURN 

END 
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SUBROUT IMF R AY 1 F ANGLE , SANGL E , XT , VT . ZT , KK ) 

THIS SUBROUTINE CALCULATES THE FCINT WHERE THE CAMERA VECTOR PIERCES 
THE SURFACE OF THE PLANET 

FANGLE - FORWARD CAMERA HALF ANGLE 
SAN3LE - SIDE CAMERA HALF ANGLE 

XT.YT.ZT - POINT WHERE THE CAMERA VECTOR PIERCES THE SURFACE 
KK - 0 IMPLIES OCES NOT PIERCE THE SURFACE-! IMPLIES DOES PIERCE 
THE SURFACE 

COMMON/BLOCK2/A.E.XI , W,0 * F , PER IOD . ANGR AT ,FC ,FTOT AL .TTCT AL 
COMMQN/BLOCK3/RD,CR,U,RS,PI,HAREF,HACOT. JDREF 

C0MMDN/fiL0CK5/IPG,IPUT,ICATE.ISTEP. ICCC . I PHOTC , I MANU . IB, IOLAP* IFR# 
♦ME, ISINGLE , 

COMMON/BLOCK8/XJ20,RES»PHI 1 , PHI 2 ,C VL AP ,C ANGF , C ANGS, FPHOTOl , FPH0T02 
♦.FOOTL , LAMBDA, L AT 1,LAT2»PLAT 
REAL LAT1 ,LAT2 ,L AMBDA,L ATL AST 
RF AL JC , JCO, JDRFF 
DIMENSION FCCTl (3,41 
KK=1 

J*I PHQTO+1 
GO TO ( 10,301 , J 

10 CONTINUE 

SF=SIN( FANGLE* DPI * CF=CCS ( F ANGLE*DR » 

SS=SIN(SANGLE*DR! * CS=COS( SANGLE*DR > 

11 CALL CCNCAR(A,E,XI,W»0»F»X»Y»Z,CX,OY*CZ,UI 
R=SORT(X*X+Y*Y+Z*ZI 

Hl=— X/R * H2*-Y/R « H3=-Z/P 

P1=H2*CZ-H3*DY $ P2=H3*DX-h 1*CZ * P3=H1*DY-H2*DX 

P=S0RT(P1*P1+P2*P2+P3*P3I 

P1=P1/P t P2*P2/P $ P3*=P3/P 

01=P2*H3-P3*H2 * Q2=P3*H1-P1*H3 $ Q3=P1*H2-P2*H1 

RH=CS*CF t RO=CS*SF * RP=SS 

R1=RH*H1+R0*01+RP*P1 

R2 = RH*H2 +P 0*02*R P*P2 

R3=RH*H3+RQ*Q3+PP*P3 

SINBC=RS/R $ COSBC=SORT (1 ,-S INBC**2 » 

IF(RH.GT.COSBC) GC TO 20 
KK=0 

GO TO 80C 

20 SR=R*RH— SORTI (R*RH! **2-R **2+R.S* *21. 

XI =X + R1*SR $ YT = Y + R 2* SR * ZT*Z+P2*SR 
GO TO 800 
30 CONTINUE 

SF=SIMFANGLE*DRl t CF= CC S < F ANCL E*OR I 

SS = SI N( ( S ANGLE + L AMBO A ) *CR ) $ CS 'COS ( ( S ANGLE+L AMBC A )*0R I 
GO TO 11 
8C0 RETURN 
END 
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SUBROUTINE STDCASE 

THIS SUBROUTINE CEFINES PRC GRAM CONSTANTS ANO STANDARD INPUT CATA 

COMMCN/PLCCK1/JDO, JD. CATE, CNUM.M.OL AST, DEL T.DELF 
COMMON/ BLCCK2/ A. E »X I t W ,0 »F » PER I CC * ANGRAT » FO » FTCT AL. TTOT AL 
CQMMQN/8LCCK3/RD»CR»U,RS,P I , F AREF » HAOOT .JDREF 
C0MM0N/BLCCK4/SX,SY,SZ,EX,EY,EZ ,CX,CY,CZ 

COMMON/ BLOCK 5/ IPG»IPUT» IC ATE, I STEP. IOCC , IPHOTO , I MANU . IB , IOLAP , IF R A 

♦ ME. I SINGLE 

COMMON/ BLOCK6/TSTOP. ESTOP , OSTOP 

COMMON/BLOCK7/HAC.HPO,AC,EC,XIC,fcC,OG.VINF,SDEC.SRA,BETA 
COMMON/BLCCKB/XJ20.RES, PHI1 ,PHI2,CVL AP,CANGF,CANGS.FPHOTOl,FPH0TC2 

♦ , FOOTL, L AMBD A. L ATI. LAT2.PL AT 
C 

REAL JD.JCn, JDREF 
REAL LATl.L AT2. LAMBDA, LATLAST 
INTEGER OSTCP.CNUM.OL AST 
DIMENSION DATEI6I»F00TL(3»4) 

C 

DATA RO. OR. PI /57. 2957795130823,. C1745 3292519943. 3. 141 59265358979/ 

C 

FT0TAL=0. 

HOT AL=0 . 

0NUM=1 
OL A S T = 1 
TOCC=Q 
IMANU=0 
IB=0 

DO 1 1=1.2 
DO 1 J-1,4 

1 FOOTHI ,J»=0. 

DO 2 J=1 *4 

2 FOOTL ( 3 , J I =1 . 

C 

TST0P=1 .E+1C 
FST0P=1.E+10 
CSTCP=10000 
U= 42E29.5 
RS= 3386. 

HAREF=149.475 
HAOOT= 350. 891962 
JDREF =24 183 22.0 
C 

XJ20=0. 00195 
RES= 10 . 

C 

RETURN 

END 
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SUBROliT I NE UPDATE! TSTEP, FSTEP) 

THIS SUBROUTINE UPDATES ALL TIME DEPENDENT PAPAMETERS SUCH AS 
TRUE ANOMALY , HOUR ANGLE , ETC. 

TSTEP - INCREASE TIME BY T STEP 
FSTEP - INCREASE TRUE ANOMALY BY FSIEP 

ANGLE (XI =AMCD(X, 360. 1+180. -SIGN! 180., XI 
C0MM0N/BL0CK1/JD0, JD , DATE ,C NUM ,F A ,0L AST . DELT , OELF 
C0 M MGN/BLCCK2/ A. E.XI , W » 0, F , PER IC D, ANGR AT ,F0 .FTOTAL , TTOT AL 
C0MMCN/BLCCK3/RD.DR.U.RS.PI ,HARE F .HADCT , JOREF 
COMMON/BLCCKA/SX.SY.SZ.EX.EY.EZ .CX .CY.CZ 

COMMON /BLOCKS/ I PG . I PUT » I C AT E . I STEF. IOCC. IPHOTO. IMANU. IB. IOLAP . IFRA 
♦ME. ISTNGLE 

COMMON/ BLCCK6/TST0P.FST OP *CSTCP 

COMMON / BLOCK 7 /H AO, HPO, AO, EC , X 10 , V>C ,00 , V I NF , SOEC . SR A, BET A 
COMMON/ BL 0CK8/XJ20.RES.PHI1.PHI2.CVLAP.CANGF.C ANGS.F PHOTO 1 .FPH0T02 
*.FnOTL,LAMROA,LATl,LAT2 .PLAT 
REAL JD.JDC.JDREF 
REAL l ATI ,L AT2.L AMBDA, LATLAST 
INTEGER OSTCP.CNUM.OL AST 
DIMENSION DATE(6).F00TL(3,4) 

JD=JD+T STEP/3600./ 24 
WDI=FLOAT( IFIXI JD) ) 

FDI= JD-WDI 

CALL JULC AL (DATE.MDI ,FDI,0 ) 

F=ANGLE( F+FSTEPl 
FTOTAL=FTCTAL+FSTEP 
TTOTAL=TTOTAL+TSTEP 
OLAST=CNUM 

ONUM=l FTOTAL +F 01 /360.+1.0 
IF ( OL AST .GT .0NUM)0LAST=0NUM 
HA=HAREF+HADCT*( JD-JDREF) 

HA= AN GL E ( F A ) 

CALL WELLS(RS .U.XJ20. A. E.XI »W.C. TSTEP I 

CALL VECTOR! JD.DECS.R AS .DECE.RAE.DECC. RAC, SX.SY.SZ.EX.EY.EZ.CX.CY, 
1CZ.4I 
RETURN 
END 
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SUBROUTINE C AL JUL ( W JO ,F JD , WND ,FC ,X) 

THIS SUBROUTINE CONVERTS A GIVEN CALENDAR DATE TO THE COPRESPC NCI KC- 
JULIAN CATE AND ALSO COMPUTES TEE NUMEER OF CAYS WHICH HAVE ELAPSEC 
SINCE JANUARY 1 AT OOCO HOURS, 195C. 

WJO - WHOLE NUMBER PART OF JUUIAA CATE 

FJO - FRACTIONAL PART OF JULIAN CATE 

WNO - WHOLE NUMBER OF CAYS SINCE JANUARY 1.0,1950 

FDA - FRACTIONAL PART OF WNC 

XI 1-61 - CALENDAR DATE ( YE AR , MONTH , CAY, HOUR , MI NUTE ,S ECONO ) 

DIMFN S ION X ( 6 ) ,AU2) 

050=243328? . 

YD=XCl)-&9. 

YL=Y0/4. 

KYL=YL 

CK=KYL 

IF( YL-CK ) 1 , 1 , 3 
1 IF(X(2)-2.)4,4,3 

3 CS=CK 
GO TO 5 

4 CS=f.K- 1 . 

5 DS=DS+365.*< YD-2.) 

DO 6 1=1,12 

6 A( I ) =1 .0 
K=X( 2) 

DO 7 I = K , 12 

7 AtlMC.O 

DS = DS + 31 . •( A (1 ) + A ( 3 ) ♦ A (5 ) + A (7 I ♦ A (8 l+A( 10 ) +4 ( 12)1 
1+30. *( AI4)+A(6)+4(9)+A( 11) )+28.*A(2) 

DS=DS+X ( 3 )-l . 

WND=DS 

FD=X (4 ) / 24 • +X I 5 ) /1440 . + X ( £ ) /864CC • 

IFCFD-. 4999999)9, 8, 8 

8 FJD=FD— .5 
W JD=1 . 

GO TO 10 

9 FJD=FD+.5 
WJD=0 . 

10 WJD=D50+WJD+WND 
RETURN 
END 


63 



uouoooouoo 


APPENDIX B 


SUBROUTINE C ARCON ( X , Y , Z , DX , OY, C Z , A » E , X I , W, 0 . F ,U I 

THIS SUBROUTINE CONVERTS CARTES IAN COORDINATES TO CONIC ELEMENTS 

X,Y.Z - COMPONENTS OF POSITION VECTOR 
DX.DY.OZ - COMPONENTS OF VELOCITY VECTOR 
A.E.XI - SEMI-MAJOR AXIS, ECCENTRICITY, INCLINATION 
W, O.F - ARGUMENT OF PERIAPSIS, LCNGITUOE OF ASCENOING NODE, TRUE 
ANOMALY 

U - GRAVITIONAL CONSTANT 

OATA RD/57.2957795130823/ 

ANGLE ( X I = 4 MOD ( X • 360 . M-18C .-SIGN ( 180. » X) 

Cl =Y*C Z— Z*DY 
C2=Z*0X-X*DZ 
C3=X*D Y— Y*DX 

H=SORT (C1*C1 ■*■02*0 2+03*03 1 
OX IR*C3/H 

XIR=ATAN2<SCRT(1.-CXIR**2) ,CXIP) 

SXIR=SIN(XIR ) 

S0=0. 

CO=l. 

IF < SXIR.EC.O.IGO TO 5 
S0=C1 /CH+SXIR) 

C0*-C2/(H*SXIR» 

5 R=S0R1( X*X+Y*Y+Z*Z) 

V=SORT I DX*DX+DY*DY+DZ*DZ ) 

A=R*U / ( 2. *U— R* V* V I 
E=S0RT(1.-H*H/(U*A»» 

P* A* C 1 ,-E*E ) 

FR=ATAN2(X*DX+Y*DY+Z*DZ»H*( P— R > / P I 

H* ANGLE ( RE* (ATAN2CI— X*SO+Y*CCI*CXIR+Z*SXIR , X*CO+Y*SO I— FR ) I 
XI=RD*XIR 

0= ANGLE (RC*ATAN2 CSO,CO> » 

F*ANGLE(RC*FR> 

RETURN 

END 
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SUBROUTINE WELL S I B.U . J2C , A ♦ E,X I ,k ,0 , OT ) 

THIS SUBROUTINE UPDATES TEE ARGUMENT OF PERIAPSIS AND THE LONGITUCE 
OF THE ASCENDING NODE 

R.U.J20 - PLANET RADIUS* GRAVITICNAL CONST ANT » 2N0 ZONAL HARMONIC 
A.E.XI - SFH I— MAJOR AXIS, ECCENTRICITY, INCLINATION 
W.O.OT - ARGUMENT OF PERIAPSIS, LONGITUDE CF ASCENDING NOOE, TIME 
INCREMENT 

REAL J20, N 

DATA PI .OP, RD/3. 14159265358979323846, .017453292519943, 
157.2957795130823/ 

N=SQRT(L/A**3) 

C=3.*MJ20*B**2 

F=A**2*(1.0-E**2>**2 

WDOT=C/F*f 1.0-15. 0/4. 0*SIN( XI*DR)**2I )*RD 

ODOT = ( I-C*COS ( X I*DR ))/ [ 2 .0 *F ) » *RC 

W=WDOT*DT+W 

C=C00T*0T*C 

RETURN 

END 
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INPUTS AND OUTPUTS 

The definitions of all input and output parameters are given in this appendix along 
with a sample input and the corresponding output. 

Definition of Input Parameters 

The following parameters are available for input into the computer program PANDG: 


Program 

symbol 

Mathematical 

symbol 

Units 

Definition 

AO 

a 

km 

Semimajor axis 

EO 

e 


Eccentricity 

XIO 

i 

deg 

Inclination 

WO 

c i ) 

deg 

Argument of periapsis 

OO 


deg 

Longitude of ascending node 

PO 

f 

deg 

Initial true anomaly 

HAO 

H a 

km 

Altitude of apoapsis 

HPO 

Hp 

km 

Altitude of periapsis 

VINF 

Vco 

km/sec 

Hyperbolic excess velocity of approach hyperbola in 
areocentric coordinate system 

SDEC 


deg 

Declination of incoming hyperbolic asymptote in areo- 
centric coordinate system 

SRA 


deg 

Right ascension of incoming hyperbolic asymptote in 
areocentric coordinate system 

BETA 

DATE 


deg 

Orientation angle of orbital plane (ref. 1) 
Six-dimensional vector defining initial calendar date. 
The order is: year, month, day, hour, minute, second 

JD 


days 

Initial Julian date 

FSTOP 

OSTOP 


deg 

Central-angle stop condition 

Orbit-number (integer) stop condition (orbit number 
changes at periapsis) 

TSTOP 


sec 

Time stop condition 

DELT 

At 

sec 

Time increment or step 

DELF 

Af 

deg 

True anomaly increment or step 

PHI1 

01 

deg 

First Sun angle denoting lighting band for photography 

PHI2 

02 

deg 

Second Sun angle 


66 



APPENDIX C 


Program 

symbol 

Mathematical 

symbol 

Units 

Definition 

FPHOTOl 

f l 

deg 

First true anomaly denoting desirable area for 




photography 

FPH0T02 

f 2 

deg 

Second true anomaly (f2 > fj) 

LATl 


deg 

First latitude denoting desirable area for photography 

LAT2 


deg 

Second latitude 

PLAT 


deg 

Latitude denoting position of a single picture 

OVLAP 



Ratio of overlap area to area of previous picture 

CANGF 

^ f 

deg 

Forward camera half-angle (see sketch (b)) 

CANGS 

^s 

deg 

Side camera half-angle (see sketch (b)) 

LAMBDA 


deg 

Out-of -plane angle for nonvertical photography (see 




sketch (b)) 

IPUT 



Input option for initial state 

1 



1 - AO, EO, XIO, WO, OO, FO 

2 - HAO, HPO, XIO, WO, OO, FO 

3 - VINF, SDEC, SRA, BETA, HAO, HPO 


IDATE 


IPG 


ISTEP 


I FRAME 


ISINGLE 


Input option for initial time 

1 - DATE (calendar date) 

2 - JD (Julian date) 

Program mode option 

1 - Groundtrack only 

2 - Photography only 

3 - Photography and groundtrack 
Step option 

1 - DELT (time increment) 

2 - DELF (true anomaly increment) 

Option denoting region of photography 

1 - Single photograph 

2 - Multiple photographs between Sun angles (PHI1, 

PHI2) 

3 - Multiple photographs between true anomalies 

(FPHOTOl, FPH0T02) 

4 - Multiple photographs between two latitudes (LATl, 

LAT2) 

Option denoting position of single photograph 

1 - PLAT (photograph taken on latitude) 

2 - Not programed 


67 



APPENDIX C 


Program 

symbol 

Mathematical 

symbol 

Units 

Definition 

IOLAP 

IPHOTO 

RS 

r s 

km 

Multiple photograph option 

0 - Successive photographs taken on step (see 

ISTEP option) 

1 - Successive photographs taken on overlap con- 

sideration (OVLAP) 

Option denoting type of photography 

0 - Vertical photography 

1 - Nonvertical photography 
Radius of planet 

U 


km^/sec^ 

Gravitational constant 

XJ20 

J 20 

None 

Second zonal harmonic of the planet 

RES 


m/km 

Photography resolution constant 

HA REF 

0o 

deg 

Hour angle at a reference Julian date 

JDREF 


day 

Reference Julian data 

HADOT 

e 

deg/day 

Time rate of change of hour angle 


Flow Diagram of Input Parameters 

A flow diagram of the input parameters has been found to be very helpful. Not all 
the parameters are defined for a specific case. Only those parameters necessary for a 
given option must be input. In addition, some numerical values have been stored in the 
program and are seldom changed. However, an option is available to change these con- 
stants when desired. Starting with the input of $CASE, the user can follow the flow dia- 
gram and choose the options or sets of input desired. An underline denotes that a numer- 
ical value must be input. For most of the parameters, the input units are given. The 
options are selected by inputing the number which is underlined. The proper input param- 
eters are chosen by proceeding through the flow diagram until the end of the case is 
reached, at which point a $ is input to terminate the case. 
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Sample Input 

The program input is loaded by using a FORTRAN IV namelist, 
input data is as follows: 

$CASE I PUT = 2, 

HAO = 30 000., 

HPO = 1000., 

XIO = 60., 

WO = 160., 

00 = 230., 

FO = 300., 

IDATE = 1, 

DATE (1) = 74., 

DATE (2) = 2., 

DATE (3) = 16., 

DATA (4) = 0., 

DATE (5) = 0., 

DATE (6) = 0., 

FSTOP = 400., 

IPG = 3, 

ISTEP = 2, DELF = 10., 

1 FRAME = 4, LAT1 = 30., LAT2 = -25, 
IOLAP = 1, OVLAP = .25, 

IPHOTO = 0, CANGF = 10, CANGS = 20.$ 


A sample set of 
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Definition of Output Parameters 
Definitions of the program output follow: 


Program 

symbol 

Mathematical 

symbol 

Units 

Definition 

JULIAN DATE 
CALENDAR DATE 

ORBIT TIME 

ORBIT NO 

JD 

day 

Current Julian date 

Current calendar date (year, month, day, hour, 
minute, second) 

Time from initial state (day, hour, minute, 
second) 

Orbit number (The orbit number is 1 initially 
and increases by 1 at each periapsis 
passage.) 

LAT 

6 

deg 

Latitude of the subsatellite point 

LONG 

a 

deg 

Longitude of the subsatellite point 

ALT 

H 

km 

Altitude of spacecraft 

VELI 

V 

km/sec 

Inertial velocity of spacecraft 

FPAI 

r 

deg 

Inertial flight-path angle 

PHI 

<p 

deg 

Sun angle (angle at subsatellite point between 
Sun vector and local vertical) 

VOHR 

V/H 

sec"f 

V over H ratio or horizontal velocity relative 
to surface divided by altitude of the 
spacecraft 

HA 

e 

deg 

Hour angle of Martian vernal equinox as mea- 
sured from the prime meridian (see 
sketch (a)) 

W 

0) 

deg 

Argument of periapsis 

O 

to 

deg 

Longitude of ascending node 

T.A. 

PHOTO OVERLAP 

FOOTPRINT 

NUMBER 

CENTER 

POINT 

f 

deg 

True anomaly of spacecraft 
Ratio of overlap area to area of previous 
footprint 

Number of footprint within present sequence 
of footprints 
Camera axis 

Points 1, 3, 5, 7 are footprint corner points 
and 2, 4, 6, 8 are at the midpoint of the sides 

LATITUDE 

6 

deg 

Latitude of respective point 

LONGITUDE 

a 

deg 

Longitude of respective point 

SUN ANGLE 

<$> 

deg 

Sun angle of respective point 

SLANT RANGE 

d 

km 

Distance between spacecraft and respective 
point (see sketch (e)) 

STATIC 

RESOLUTION 


m 

Static resolution is slant range multiplied by a 
constant (RES) 
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Sample Output 

The output which resulted from the sample input described previously is presented 
subsequently. These output data are also presented graphically in figure 1. Because of 
the lack of space, only the first part of the computer printout follows. 
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